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ABSTRACT 

We present deep optical and near-infrared (NIR) UBVRIHKs imaging data for 24 
blue compact galaxies (BCGs). The individual exposure times are on average ~ 40 
minutes in the optical (B) and ~ 90 minutes in the NIR, but on occasion up to ~ 5 
hours for a single target and filter, observed on 2.5, 3.5, 8.2— m telescopes. The sam- 
ple contains luminous dwarf and intermediate-mass BCGs which are predominantly 
metal-poor, although a few have near-solar metallicities. We have analyzed isophotal 
and elliptical integration surface brightness and color profiles, extremely deep (ps ^29 
mag arcsec -2 ) contour maps and RGB images for each galaxy in the sample, and pro- 
vide a morphological classification where such is missing. We have measured the total 
galaxy colors, the colors of the underlying host galaxy, and the colors of the burst, and 
compare these to the predictions of new state-of-the-art spectral evolutionary mod- 
els (SEMs) both with and without contribution by nebular emission. Separating the 
burst from the underlying host we find that regardless of the total luminosity the host 
galaxy has the properties of a low surface brightness (LSB) dwarf with Mb ^ —18. 
For a number of galaxies we discover a distinct LSB component dominant around and 
beyond the Holmberg radius. For the specific case of ESO400-43A&B we detect an 
optical bridge between the two companion galaxies at the fly ~ 28th mag arcsec -2 
isophotal level. Synthetic disk tests are performed to verify that we can trace such 
faint components with negligible errors down to \ib = 28 mag arcsec -2 and \ik = 23 
mag arcsec -2 . By examining the structural parameters (central surface brightness /j,q 
and scale length h r ) derived from two radial ranges typically assumed to be dominated 
by the underlying host galaxy, we demonstrate the importance of sampling the host 
well away from the effects of the burst. We find that /zo and h r of the BCGs host 
deviate from those of dwarf ellipticals (dE) and dwarf irregulars (dl) solely due to a 
strong burst contribution to the surface brightness profile almost down to the Holm- 
berg radius. Structural parameters obtained from a fainter region, /is = 26-28 mag 
arcsec -2 , are consistent with those of true LSB galaxies for the starbursting BCGs in 
our sample, and with dEs and dls for the BCGs with less vigorous star formation. 
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1 INTRODUCTION 

Blue compact galaxies (BCGs) are gas-rich star-forming 
low redshift galaxies with low metallicities (for a review 
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see iKunth fc Ostlinl |200Ch . Their star formation rates 
(SFRs) can be as high as ~ 20 M© yr -1 (e.g. Haro 11, 
lHaves et alj 120071 ) and their H I gas masse s range from 
M(HI) ~ 10 7 to a few t ens of 10 8 Mq (e.g. iTavlor et all 
1 19951 ; iThuan et all 1 19991 ; ISmoker et at] |2000D . BCGs are 
said to be starbursting when they will exhaust their gas 
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supply in less than a Hubble time. The metallicities mea- 
sured in these galax i es are low (e.g. [ Masegosa et al.l 1 19941 ; 
iMarconi et~aT1 1 19941; flzotov fc Thuani 1 19991 )," to extremely 
low l lzotov et al.l 20051 ), once leading rese archers to propose 
that BCGs are truly young galaxies (|Searle fc Sargentl 
1 19721 ), experiencing their first episode of star formation in 
a primordial gas environment. There is mounting evidence 
to the contrary, however, and "first-burst" claims have 
been cha llenged using broadband colors and morphol- 
ogy (e.g. IPapaderos et all Il996l; iTelles fc Terlevichl 1 19971 ; 



ICairos et al.l l2001al lbl: iBergvall fc Ostlinl |2002|), a bsorption 
features in the spectra (e.g. IWestera et al.ll2004l ). through 
the di rect detection of gl obular clusters associated with the 
BCG (|Qstlin et al.lll998l ). and for near-by targets through 
the detection of luminous red giants and asymptotic giant 
branch (AGB) s t ars in color-m agnitude diagrams (e.g. 
lAloisi et ail 1 19991 ; lOstlml l2000al ). In most cases the rel- 
ative strength of the starburst can be so high that it 
completely dominates the light output of the galaxy, 
an obstacle which has been countered by deeper optical 
imaging dat a and observations i n the near infra-re d (NIR) 
regime (e.g. iNoeske et all 120031 ; ICairos et all 120031 ) and an 
older population has been revealed. In other cases, e.g. 
SBS0335-052 E&W, the burst can be spatially dominating 
over the galaxy even for extremely deep imaging, making 
it difficult to directly detect a possible old population. 
Hence, the youth hypo t hesis persists for th ese targets (e.g. 
iThuan fc Izotovl Il997l: llzotov et al.l 120051 ) but has been 
challenged (|6stlin fc Kunthll200ll ) 



Nevertheless, the emerging consensus seems to be that 
BCGs are galaxies with an underlying old stellar population, 
often referred to as 'host', on which the current starburst 
is superposed. The host population is distinctly different 
in te r ms of optical and NIR colors (e . g. Bergvall fc Ostlinl 
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2002]; INoeske et al.ll2003l; [Cairos et al 
2005; Caon et al.ll2005l ; I Amorin et al.ll2009l ). and appears to 
be smoothly distributed in a phy sical structure rem i niscent 
of a disk with regular isophotes (lLoose et al.lll983l ; iThuanl 
1 19831 ; IPapaderos et all Il996l ; ISmoker et al.lll999l ). The dis- 
covery of this host population invariably led to a relaxation 
in the definition of this group of galaxies. The "dwarf" 
qualifier in the original naming convention of "blue compact 
dwarf" is often dropped by some authors to the benefit of 
the more inclusive "galaxy" term. There is a large fraction 
of BCGs which are too massive and bright to be dwarfs, 
with Mb as high as —20, but they too exhibit starbursts in 
a gas-rich, metal-poor environment, and a similar morpho- 
logical structure as the dwarfs - an extended, often regular, 
disk-like structure hosts knots of star formation (SF), 
either centrally located or dispersed in seemingly random 
fashion around the gal axy. The name i s still misleading 
since, as pointed out in lGil de Paz et all (|2003l ). some blue 
compact galaxies are neither blue, nor compac t . Incl uding 
the host, the requirements of lThuan fc Martini (|l98ll ) that 
a BCG must have an apparent optical size 7-25 < 1 kpc 
in diameter can now only be met for some very compact 
sources, like HL293B. Further, the "bl ue" condition of 
B - V < 0.3 m (|Bergvall fc 01ofsson|[l98rj ) can now only be 
applied to the starbursting region since inclusion of the host 
shifts the total B — V color redwards of this (somewhat 



arbitrary) limit. 

Despite such diversity all BCGs have at least one thing in 
common. They are a subgroup of emission line galaxies, 
well separated from Seyferts and quasars, with spectra 
reminiscent of H II regions, only with a strong stellar con- 
tinuum from the composite populations of th e galaxy (e.g. 
Kunth fc Sargentl 1 19831 ; iMasegosa et al.l Il994l ; ISalzer et all 
1989a! ) . Because of their H II region-like spectrum they are 
often referred to as H II galaxies. The SFR can drastically 
vary from BCG to BCG, with some extreme outliers like 
Haro 11 having a SFR > 18 Mp, yr -1 , which it can 
only sustain for < 50 Myr (|Bergvall et all |2000| ). While 
not equally vigorous, the SF activity in most BCGs must 
similarly proceed in short-lived bursts, lasting from a few 
tens of Myr to 1 Gyr. The duration of the typical burst has 
been constrained both from observed H I masses, which 
are generally too low to sustain the current SFRs for a 
Hubble time, and from measured metallicities, which are 
also too low for the interstellar medium to have undergone 
enrichment by constant star formation. The latter argument 
can be partially circumvented if one allows for the outflow 
and loss of metals to the intergalactic medium due to 
supernovae winds, or for the inflow of pristine gas, which 
would dilute the metallicity. The argument of too low 
metallicity cannot be completely discarded, however, as 
there are indications that gas-rich galaxies ar e more likely 
to be closed systems (|van Zee fc Havnesll2006l ). 

Some BCGs are clearly in various stages of merging, or 
undergoing tidal interaction. They can have very irregular 
morphologies, with complicated, almost chaotic kinematics, 
making it difficult to quantify the masses and ages of 
the current burst. Other BCGs show an otherwise reg- 
ular elliptical disk with a central nuclear starburst, or 
sporadic star-for ming knots offset f rom the center, a.k.a. 
H II hotspots l|Salzer et al.l Il989al ), and it is not clear 
whether galaxies with such morphologies also are subject 
to the same SF trigger. H I observations have shown that 
while the majority of BCGs are associated with massive 
H I clouds, most likely fueling the ongoing star formation, 
they are generally isolated in the sense that they lack 
mass iv e companion galaxie s dCampos-Aguilar fc Molesl 
199ll; ICampos-Aguilar et all Il993l ; iTavlor et al.l 1 19951 : 
Thuan et all 1 19991 ). However, H I companions without 
a co nfirmed optic al counterpart do appear to be com- 
mon l|Tavloijll997l ). indicating that perhaps the companion 
galaxy is a faint gas-rich dwarf. This makes mergers the 
leading candidates for SF triggers since tidal interactions 
are unlikely to be the primary SF triggers in dwarf galax- 
ies l|Brosch et all 12004). For a starburst to take place the 
merger would not have to be major but it would have to 
be wet, meaning at least one of the merging galaxies would 
have to be gas-rich. 

Ideally, one would like to study the SF as it occurs for 
the first time in the protogalaxies in the distant Universe. 
It is thought that observing such pristine SF activity 
would allow us to determine what the necessary trigger 
conditions for the onset of star formation are. Presently, 
we are unable to directly study the stellar populations 
of high-redshift galaxies in any great detail since we are 
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constrained by cosmological effects and small intrinsic 
sizes of the targets, and the comparatively poor resolution 
of our currently available instrumentation. Because of 
their low metallicities, high SFRs, and their lack of spiral 
density waves, BCGs are thought to be reminiscent of 
the young galaxies in the universe, in which the star 
formation must have proceeded under similar metal- 
poor gas-rich conditions. Being much closer, BCGs are 
available for detailed photometric and spectroscopic studies. 

Considering the degree of complication inherent to the anal- 
ysis of such a heterogeneous group of galaxies, one needs to 
examine the properties of the BCGs from several angles: the 
behavior of the spectrum in the optical and NIR regimes 
on a large scale (broadband imaging), the individual line 
emissions (typically H a and [Oil I]), and the kinematics 
of the gas and the stars. We have undertaken a study 
of BCGs that combines optical and NIR broadband a nd 
narrowband imaging (this work, Mi cheva et al.l l2012al lbh . 
Fabry-Perot interferometry (Ostlin et al. 2012 in prep., 
Marquart et al. 2012 in prep), and grism spectroscopy of 
two distinct types of star-forming galaxies - high and low 
luminosity BCGs. In this paper, which is the first in a series, 
we present broadband UBVRIHKs data for a sample 
of 24 luminous galaxies classified as BCGs. This sample 
was hand-picked to contain interesting and representative 
cases, and is biased towards relatively luminous (median 
M b y —18 mag) gala xies. Some of the target s come from 
the Ostlin et al.l |200ll ) sample, others from the lSalzer et al.l 
(|l989al ) UM survey of H II galaxies. The latter were 
selected based on their emission line spectra to be actively 
star-forming if not always strongly starbursting, and to 
have metallicities and H I masses available in the literature. 
Additionally, a number of famous and strongly starbursting 
BCGs (e.g. IIZw40, MK930, Haroll) were included in the 
project because they lacked equally deep observations and 
yet presented interesting targets to have in our analysis . 
Haro 11 was already presented in iMicheva et"aH l|2010l ). 
so it is not presented here but will be in cluded in our 
future analysis. In a consecutive paper (jMicheva et al.l 
l2012al ) we present UBVRIHKs broadband imaging data 
for a volume-limited sample of 21 low luminosity emission 
line galaxies. These constitute a well defined complete 
sample of star forming galaxies, representative of the galaxy 
population in the local Universe. In total we have 46 BCG 
and BCG-like targets in 7 broadband and 2 narrowband 
filters. In the optical, imaging observations of depth com- 
parable to ou rs can be found in the li terature for a number 
of galaxies jlPapaderos et al. _ 19961; Doublier et alj 1 19991 ; 
ICairos et all l2001al ; iBergvall fc Ostlinl |2002l ). In "the NIR, 
however, there is no sample of BCGs wit h as deep individual 
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exposure times as our s (compare e.|^ _ 
ICairos et al1l2003l ). In IMicheva et alj (2012b) we combine 
both samples to juxtapose the high and low luminosity 
BCGs and analyse any differences in their host populations. 
In the latter paper we also present the narrowband data 
from both samples as well as derived estimates of the ages, 
masses and star formation rates through SED fitting. We 
have assumed Ho = 73 km -1 s _1 Mpc -1 throughout this 
paper. 



servations, reduction and calibration of the data. §|3j §|4] and 
§ [5] deal with the generation of surface brightness profiles, 
contour plots and three-color RGB images, respectively. § [5] 
and §[7]explain our approach to obtaining total luminosities 
and total colors for the galaxies. §|8]illustrates how we obtain 
structural parameters for population separation. Character- 
istics of individual galaxies and comparison of integrated 
colors to stellar evolutionary models can be found in § 1111 
Finally, discussion and concluding remarks are in § and 

sua 



2 OBSERVATIONS 

The data consist of optical and near infra-red (NIR) broad- 
band imaging in the U BVRIHK filters. They were obtained 
during the period 2001-2007 with ALFOSC (Nordic Optical 
Telescope, NOT), EMMI (European Southern Observatory 
New Technology Telescope, ESO NTT), and FORS (ESO 
Very Large Telescope, VLT) in the optical, and with NOT- 
CAM (NOT) and SOFI (ESO NTT) in the NIR. A log of 
observations is shown in Tableland the individual exposure 
times for each filter in Table [1] Since both t he sample pre- 
sented here and the volume limited sample of lMicheva et al.l 
<|2012al ) were reduced with the same pipelines, the descrip- 
tion of the reductions and calibration that follows here ap- 
plies to both. We will give the details in t his paper and will 
not repeat them in IMicheva et all (|2012aT ). 



The layout of this paper is as follows: §[2] introduces the ob- 
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Table 2. Log of the observations. Heliocentric redshift and cosmology-corrected luminosity distances from NED. 



Galaxy Ra Dec (J2000) Redshift D [Mpc] Year Instrument Filters 



ES0185- 


13 


19h45m 


00.5a 


0. 


018633 


76 


3 


2001 


EM MI 


Bb#605,V#606,I#610 






-54dl5r 


n03s 












SOFI 


Ks#13 






















U BESS R BESS 




















SOFI 


H#12 


ES0249- 


31 


03h55m 


46.1s 


0. 


002829 


10 


8 


2001 


EMMI 


Bb#605,V#606,R#608,I#610 






-42d22r 


n05s 












SOFI 


Ks#13 






















H#12 


ES033S- 


104 


19h27m 


58. 2s 


0. 


009453 


37.4 


2001 


EMMI 


R#608 






-41d34r 


n32a 












EMMI 


Bb#605,V#606 




















EMMI 




















2005 


SOFI 


Ks#13 


















2006 


SOFI 


Ks#13,H#12 




















E A I M I 


\ 1 j 1 1 1 J .1 ^fp o ± u 






-35d29r 


n08s 










2006 


SOFI 


KS#13,H#12 


ES0421- 


02 


04h29m 


40.1s 


0. 


003149 


12 




2001 


EMMI 


Bb#605,V#606,R#608,I#610 






-27d24r 


ii31s 












SOFI 


Ks#13 


















2002 


FORS1 


U_BESS 




















SOFI 


H#12 


ES0462- 


IG020 


20h26m 


56.8s 


0. 


019782 






2001 


EMMI 


Bb#605,V#606 






-29d07r 


nOls 












SOFI 


Ks#13 


















2002 


FORS1 


U.BESS 




















SOFI 


H#12 


HE2-10 




08h36m 


15.1s 


0. 


002912 


15 


8 


2001 


EMMI 


Bb#605,V#606 






-26d24r 


n34s 










2002 


FORS1 


U_BESS 




















SOFI 


H#12,Ks#13 


HL293B 




22h30m 


36.798a 





005170 


16 


3 


2001 


ALFOSC-FASU 


B#74,V#75,R#76,I#12 






-00d06r 


q36.99s 










2002 


ALFOSC-FASU 


U#7 


















2003 


NOTCAM 


Ks#207,H#204 


IIZW40 




05h55m 


42.6s 


0. 


002632 


11. 


8 


2001 


ALFOSC-FASU 


B#74,I#12 






+03d23 


m32s 










2002 


ALFOSC-FASU 


U#7 




















NOTCAM 


Ks#207 




















SOFI 


Ks#13 


















2003 


ALFOSC-FASU 


V#75 




















N OTC AM 


H # 204 


MK600 




02h51m 


04.6s 


0. 


003362 


10 


9 


2001 


ALFOSC-FASU 


B#74,I#12 






+04d27 


ml4s 










2002 


ALFOSC-FASU 


U#7,V#75,R#76 




















N OTC AM 




















2003 


NOTCAM 


H#204 


MK900 




21h29m 


59.6s 


0. 


003843 


. ! 




2001 


ALFOSC-FASU 


B#74,V#75,R#76,I#12 






_j_Q2d24 














ALFOSC-FASU 




















2003 


NOTCAM 


Ks#207,H#204 


MK930 








0. 


018296 




' 




ALFOSC-FASU 


T3-U.-7A \r-U--7K TJIIO 

L3 =fp i <-± , V =fp i O t^-ii 1 J- ^ 






+ 28d56 


m50s 










2002 


ALFOSC-FASU 


U#7,R#76 


















2003 


NOTCAM 


Ks#207,H#204 


MK996 








0. 


005410 


18 


2 




ALFOSC-FASU 


I#12 






-06dl9r 


q36s 










2002 


ALFOSC-FASU 


U#7 




















SOFI 


Ks#13 


















2003 


ALFOSC-FASU 


B#74,V#75,R#76 




















NOTCAM 


H#204 


SBS0335 


-052E&W 


03h37m 


44.0a 


0. 


013486 


54 




2001 


ALFOSC-FASU 


B#74,I#12 






-05d02r 


n40s 










2002 


SOFI 


Ks#13 


















2003 


NOTCAM 


KS#207,H#204 


TOL0341 


03h42m 


49.4s 


0. 


014827 


60 


6 


2001 


EMMI 


Bb#605,V#606,I#610 






-40d35r 


n56s 










2002 


FORS1 


U_BESS 




















SOFI 


H#12,Ks#13 


TOL1457-2627&77 


15h00m 


27.9s 


0.016812 


72 


.8 


2004 


ALFOSC-FASU 


B#74,V#75 






-26d27r 


n02s 










2005 


SOFI 


Ks#13 


UM133 




01h44m 


41.3s 


0. 


005414 


18. 


3 


2001 


ALFOSC-FASU 


B#74,I#12 






+04d53 


m26s 










2002 


ALFOSC-FASU 


U#7,V#75,R#76 


















2003 


NOTCAM 


Ks#207,H#204 


UM160 




23h24m 


19.8s 


0. 


008006 


28 




2001 


ALFOSC-FASU 


B#74,V#75,I#12 






-00d07r 


nOls 










2002 


ALFOSC-FASU 


U#7,R#76 




















SOFI 


H#12,Ks#13 


UM238 




00h24m 


42.3a 


0. 


014230 


54 


2 


2001 


ALFOSC-FASU 


B#74,I#12 






+01d44 


m02s 










2002 


ALFOSC-FASU 


U#7,V#75,R#76 




















SOFI 


Ks#13 


















2003 


NOTCAM 


H#204 


UM417 




02hl9m 


30.2a 


0. 


009000 


33 


8 


2001 


ALFOSC-FASU 


B#74,I#12 






-00d59r 


nils 










2002 


ALFOSC-FASU 


U#7,V#75,R#76 



NOTCAM Ks#207 
2003 NOTCAM H#204 



UM448 


Ilh42ml2.4a 


0.018559 


74.6 


2004 


ALFOSC-FASU 


B#74,V#75 




4-00d20m03s 






2005 


ALFOSC-FASU 


U#7 












EMMI 


R#608,I#610 










2007 


NOTCAM 


Ks#207 


UM619 


13h52m44.7s 


0.015494 


63.1 


2004 


ALFOSC-FASU 


B#74,V#75 




+00d07m53a 






2005 


ALFOSC-FASU 


U#7 












EMMI 


R#608,I#610 












NOTCAM 


H#204 












SOFI 


Ks#13 










2006 


SOFI 


Ks#13 










2007 


NOTCAM 


Ks#207 



© 2011 RAS, MNRAS OOO.ITU621 



Deep surface photometry on 24 BCGs 



Radius [Kpc] 

4 6 13 



K 




Radius farcsecl 

Figure 1. ES0185-13. Left panel: Surface brightness and color radial profiles for elliptical (open circles) and 
isophotal (red circles) integration. Upper right panel: contour plot based on the B band. Isophotes fainter than 
22.5, 25.5 are iteratively smoothed with a boxcar median filter of size 5, 15 pixels respectively. Lower right panel: A 
true color RGB compos i te im age using the U, B, I filters. Each channel has been corrected for Galactic extinction 
following ISchlegel et al.l dl998h and co nverted to the AB photometric system. The RGB composite was created by 
implementing the lLupton et aL | d2004h algorithm. 
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Table 1. Total integration times for the sample. All 
times are given in minutes and converted to the frame- 
work of a 2.56 meter telescope where needed. The 
values are for observations in a single filter, e.g. only 
SOFI Ks, and not SOFI Ks + NOTCAM Ks. 
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Figure 1 — continued ES0249-31. Left panel: Surface brightness and color radial profiles for elliptical (open circles) 
and isophotal (red circles) integration. Upper right panel: contour plot based on the B band. Isophotes fainter than 
22.0, 24.5, 25.5 are iteratively smoothed with a boxcar median filter of size 5, 15, 25 pixels respectively. Lower right 
panel: A true color RGB composite ima ge using the B, V, I filters. Each channel has been corrected for Galactic 
extinction following ISchlegel et al.l dl998f) and co nverted to the AB photometric system. The RGB composite was 
created by implementing the lLupton et al.l ll2004r ) algorithm. 
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Figure 1 — continued ESO338-04. Left panel: Surface brightness and color radial profiles for elliptical (open circles) 
and isophotal (red circles) integration. Upper right panel: contour plot based on the B band. Isophotes fainter than 
22.5, 25.5 are iteratively smoothed with a boxcar median filter of size 5, 15 pixels respectively. Lower right panel: A 
true color RGB compos i te im age using the B,V,I filters. Each channel has been corrected for Galactic extinction 
following ISchlegel et al] lll998h and co nverted to the AB photometric system. The RGB composite was created by 
implementing the lLupton et al" | d2004h algorithm. 
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Figure 1 — continued ESO400-43. Left panel: Surface brightness and color radial profiles for elliptical (open circles) 
and isophotal (red circles) integration. Upper right panel: contour plot based on the V band. Isophotes fainter than 
23.5, 25.5 are iteratively smoothed with a boxcar median filter of size 5, 15 pixels respectively. Lower right panel: A 
true color RGB compos i te im age using the V, I, H filters. Each channel has been corrected for Galactic extinction 
following ISchlegel et al.l dl998h and co nverted to the AB photometric system. The RGB composite was created by 
implementing the lLupton et ah | d2004h algorithm. 
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Figure 1 — continued ESO421-02. Left panel: Surface brightness and color radial profiles for elliptical (open circles) 
and isophotal (red circles) integration. Upper right panel: contour plot based on the B band. Isophotes fainter than 
23.5, 25.5 are iteratively smoothed with a boxcar median filter of size 5, 15 pixels respectively. Lower right panel: A 
true color RGB compos i te im age using the U, B, I filters. Each channel has been corrected for Galactic extinction 
following ISchlegel et al] dl998h and co nverted to the AB photometric system. The RGB composite was created by 
implementing the lLupton et ah | d2004h algorithm. 
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Figure 1 — continued ESO462-20. Left panel: Surface brightness and color radial profiles for elliptical (open circles) 
and isophotal (red circles) integration. Upper right panel: contour plot based on the B band. Isophotes fainter than 
22.5 are smoothed with a boxcar median filter of size 5 pixels. Lower right panel: A true co lor RGB composite i mage 
using the U, B, V filters. Each channel has been corrected for Galactic extinction following lSchle gel et al. l jl998l') and 
converted to the AB photometric system. The RGB composite was created by implementing the lLuptmietaTl |2004l ) 
algorithm. 
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Figure 1 — continued HE2-10. Left panel: Surface brightness and color radial profiles for elliptical (open circles) 
and isophotal (red circles) integration. Upper right panel: contour plot based on the B band. Isophotes fainter than 
22.5, 24.5, 25.5 are iteratively smoothed with a boxcar median filter of size 5, 15, 25 pixels respectively. Lower right 
panel: A true color RGB composite imag e using the U, B, V filters. Each channel has been corrected for Galactic 
extinction following Schl egel et al.| | l998fl and co nverted to the AB photometric system. The RGB composite was 
created by implementing the lLupton et al. algorithm. 
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Figure 1 - continued HL293B. Left panel: Surface brightness and color radial profiles for elliptical (open circles) 
and isophotal (red circles) integration. Upper right panel: contour plot based on the B band. Isophotes fainter than 
23.5, 25.5 are iteratively smoothed with a boxcar median filter of size 5, 15 pixels respectively. Lower right panel: A 
true color RGB compos i te im age using the U, B, I filters. Each channel has been corrected for Galactic extinction 
following ISchlegel et al.l dl998h and co nverted to the AB photometric system. The RGB composite was created by 
implementing the lLupton et al" | d2004h algorithm. 
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Figure 1 — continued IIZw40. Left panel: Surface brightness and color radial profiles for elliptical (open circles) 
and isophotal (red circles) integration. Upper right panel: contour plot based on the B band. Isophotes fainter than 
22.5, 25.0 are iteratively smoothed with a boxcar median filter of size 5, 15 pixels respectively. Lower right panel: A 
true colo r RGB compos i te im age using the B, V, I filters. Each channel has been corrected for Galactic extinction 
following Schlegel et al.| lj l998h and co nverted to the AB photometric system. The RGB composite was created by 
implementing the lLupton et ah | d2004h algorithm. 
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Figure 1 — continued MK600. Left panel: Surface brightness and color radial profiles for elliptical (open circles) 
and isophotal (red circles) integration. Upper right panel: contour plot based on the B band. Isophotes fainter than 
23.0, 25.0 are iteratively smoothed with a boxcar median filter of size 5, 15 pixels respectively. Lower right panel: A 
true color RGB compos i te im age using the U, B, I filters. Each channel has been corrected for Galactic extinction 
following ISchlegel et al.l dl998h and co nverted to the AB photometric system. The RGB composite was created by 
implementing the lLupton et ah | d2004h algorithm. 
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Figure 1 — continued MK900. Left panel: Surface brightness and color radial profiles for elliptical (open circles) 
and isophotal (red circles) integration. Upper right panel: contour plot based on the B band. Isophotes fainter than 
22.5, 25.5 are iteratively smoothed with a boxcar median filter of size 5, f5 pixels respectively. Lower right panel: A 
true color RGB compos i te im age using the U, B, I filters. Each channel has been corrected for Galactic extinction 
following ISchlegel et al] dl998h and co nverted to the AB photometric system. The RGB composite was created by 
implementing the lLupton et al" | d2004h algorithm. 
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Figure 1 — continued MK930. Left panel: Surface brightness and color radial profiles for elliptical (open circles) 
and isophotal (red circles) integration. Upper right panel: contour plot based on the B band. Isophotes fainter than 
24.0, 26.0 are iteratively smoothed with a boxcar median filter of size 5, f5 pixels respectively. Lower right panel: A 
true color RGB compos i te im age using the U, B, I filters. Each channel has been corrected for Galactic extinction 
following ISchlegel et al.l dl998h and co nverted to the AB photometric system. The RGB composite was created by 
implementing the lLupton et ah | d2004h algorithm. 
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Figure 1 — continued MK996. Left panel: Surface brightness and color radial profiles for elliptical (open circles) 
and isophotal (red circles) integration. Upper right panel: contour plot based on the B band. Isophotes fainter than 
23.0, 25.5 are iteratively smoothed with a boxcar median filter of size 5, 15 pixels respectively. Lower right panel: A 
true color RGB compos i te im age using the U, B, I filters. Each channel has been corrected for Galactic extinction 
following ISchlegel et al] dl998h and co nverted to the AB photometric system. The RGB composite was created by 
implementing the lLupton et ah | d2004h algorithm. 
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Figure 1 — continued SBS0335-052 EAST. Left panel: Surface brightness and color radial profiles for elliptical 
(open circles) and isophotal (red circles) integration. Upper right panel: contour plot based on the B band. Isophotes 
fainter than 23.5, 26.0 are iteratively smoothed with a boxcar median filter of size 5, 15 pixels respectively. Lower 
right panel: A true c olor RGB composite i mage using the B, I, K filters. Each channel has been corrected for Galactic 
extinction following ISchlegel et al.l jl99ST) and co nverted to the AB photometric system. The RGB composite was 
created by implementing the lLupton et al.l (120041 ') algorithm. 
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Figure 1 - continued SBS0335-052 WEST. Le/i panel: Surface brightness and color radial profiles for elliptical 
(open circles) and isophotal (red circles) integration. Upper right panel: contour plot based on the B band. Isophotes 
fainter than 23.5, 26.0 are iteratively smoothed with a boxcar median filter of size 5, 15 pixels respectively. Lower 
right panel: A true c olor RGB composite i mage using the B, I, K filters. Each channel has been corrected for Galactic 
extinction following Schlegel et al.| | l998fl and co nverted to the AB photometric system. The RGB composite was 



created by implementing the lLupton et al. algorithm. 
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Figure 1 — continued TOL0341-407. Left panel: Surface brightness and color radial profiles for elliptical (open 
circles) and isophotal (red circles) integration. Upper right panel: contour plot based on the B band. Isophotes 
fainter than 23.5, 25.5 are iteratively smoothed with a boxcar median filter of size 5, 15 pixels respectively. Lower 
right panel: A true c olor RGB composite i mage using the U, B, I filters. Each channel has been corrected for Galactic 
extinction following ISchlegel et al.l dl998f) and co nverted to the AB photometric system. The RGB composite was 
created by implementing the lLupton et al.l (120041 ') algorithm. 
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Figure 1 — continued TOL1457— 262 7. Left panel: Surface brightness and color radial profiles for elliptical (open 
circles) and isophotal (red circles) integration. Upper right panel: contour plot based on the B band. Isophotes fainter 
than 24.0, 25.5 are iteratively smoothed with a boxcar median filter of size 5, 15 pixels respectively. Lower right 
panel: A true color RGB composite imag e using the B,V,K filters. Each channel has been corrected for Galactic 
extinction following Schlegel et al.| | l998fl and co nverted to the AB photometric system. The RGB composite was 
created by implementing the lLupton et al. algorithm. 
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Figure 1 — continued TOL1457— 262 II. Left panel: Surface brightness and color radial profiles for elliptical (open 
circles) and isophotal (red circles) integration. Upper right panel: contour plot based on the B band. Isophotes fainter 
than 24.0, 25.5 are iteratively smoothed with a boxcar median filter of size 5, 15 pixels respectively. Lower right 
panel: A true color RGB composite imag e using the B,V,K filters. Each channel has been corrected for Galactic 
extinction following Schlegel et al.| | l998fl and co nverted to the AB photometric system. The RGB composite was 
created by implementing the lLupton et al. algorithm. 
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Figure 1 — continued UM133. Left panel: Surface brightness and color radial profiles for elliptical (open circles) 
and isophotal (red circles) integration. Upper right panel: contour plot based on the B band. Isophotes fainter than 
23.0, 25.5 are iteratively smoothed with a boxcar median filter of size 5, 15 pixels respectively. Lower right panel: A 
true color RGB compos i te im age using the U, B, I filters. Each channel has been corrected for Galactic extinction 
following ISchlegel et al] dl998h and co nverted to the AB photometric system. The RGB composite was created by 
implementing the lLupton et ah | d2004h algorithm. 
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Figure 1 — continued UM160. Left panel: Surface brightness and color radial profiles for elliptical (open circles) 
and isophotal (red circles) integration. Upper right panel: contour plot based on the B band. Isophotes fainter than 
23.5 are smoothed with a boxcar median filter of size 5 pixels. Lower right panel: A true c olor RGB composite i mage 
using the U, B, I niters. Each channel has been corrected for Galactic extinction followine lSchle gel et al. l dl99St) and 
converted to the AB photometric system. The RGB composite was created by implementing the lLuptonetaTI l|2004f ) 
algorithm. 



© 2011 RAS, MNRAS OOO.ITH621 



Micheva et al. 



Radius [kpc] 

3 6 -ft 




5 10 15 20 25 30 35 40 
Radius larcsecl 



Figure 1 — continued UM238. Left panel: Surface brightness and color radial profiles for elliptical (open circles) 
and isophotal (red circles) integration. Upper right panel: contour plot based on the B band. Isophotes fainter than 
23.5, 25.5 are iteratively smoothed with a boxcar median filter of size 5, 15 pixels respectively. Lower right panel: A 
true color RGB compos i te im age using the U, B, I filters. Each channel has been corrected for Galactic extinction 
following ISchlegel et al] dl998h and co nverted to the AB photometric system. The RGB composite was created by 
implementing the lLupton et ah | d2004h algorithm. 
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Figure 1 — continued UM417. Left panel: Surface brightness and color radial profiles for elliptical (open circles) 
and isophotal (red circles) integration. Upper right panel: contour plot based on the B band. Isophotes fainter than 
23.5, 26.0 are iteratively smoothed with a boxcar median filter of size 5, 15 pixels respectively. Lower right panel: A 
true color RGB compos i te im age using the U, B, I filters. Each channel has been corrected for Galactic extinction 
following ISchlegel et al.l dl998h and co nverted to the AB photometric system. The RGB composite was created by 
implementing the lLupton et ah | d2004h algorithm. 
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Figure 1 — continued UM448. Left panel: Surface brightness and color radial profiles for elliptical (open circles) 
and isophotal (red circles) integration. Upper right panel: contour plot based on the B band. Isophotes fainter than 
23.9, 25.5 are iteratively smoothed with a boxcar median filter of size 5, 15 pixels respectively. Lower right panel: A 
true color RGB compos i te im age using the B,V,I filters. Each channel has been corrected for Galactic extinction 
following ISchlegel et al] dl998h and co nverted to the AB photometric system. The RGB composite was created by 
implementing the lLupton et ah | d2004h algorithm. 
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Figure 1 — continued UM619. Left panel: Surface brightness and color radial profiles for elliptical (open circles) 
and isophotal (red circles) integration. Upper right panel: contour plot based on the B band. Isophotes fainter than 
23.5, 25.5 are iteratively smoothed with a boxcar median filter of size 5, 15 pixels respectively. Lower right panel: A 
true color RGB compos i te im age using the U, B, I filters. Each channel has been corrected for Galactic extinction 
following ISchlegel et al.l [jl998h and co nverted to the AB photometric system. The RGB composite was created by 
implementing the lLupton et al" | d2004h algorithm. 
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2.1 Data reductions 

In order to ensure the data consistently undergo the 
same processing, we have reduced all raw frames, except 
BVR data for ESO338-04, with two pipelines specifically 
designed for this project - one for the optical and one for the 
NIR. This consistency ensures us that any difference in the 
results can not be attributed to the reduction process itself 
and instead reflects a true variation in properties between 
different galaxies. Both pipeli nes were already described in 
detail in lMicheva et al.l ((2010) , so here we will only present 
additional points of interest relevant to the reduction. 

The reductions were carried out on subsamples grouped by 
night and instrument. For every night and all instruments 
an individual sky background was fit to each reduced (i.e. 
bias-subtracted and flatnelded) frame before combining 
them to a final stack, and at no point in the reduction 
process was a "master sky" ever used. After stacking, a 
final sky fitting and subtraction was performed. Calibration 
in both wavelength regimes is automatic, using Landolt 
standard stars in the optical and Two Micron All Sky 
Survey (2MASS) stars in the NIR. 

ALFOSC Some nights on this instrument were lacking 
usable standard stars. However, for all such observations 
short 3 minute exposures were taken during photometric 
nights, or the same target was observed on several occasions. 
It was therefore always possible to calibrate problematic 
nights against photometric observations of the same object. 

For 2004 data there was a severe charge bleed in the frames 
of UM439, which the observers attempted to fix by varying 
the rotator field angle between consecutive frames of that 
target. This varying orientation of the images proved to be 
too much for the IRAF XYXYMATCH procedure, which 
the pipeline uses on SExtractor sources in order to align 
the frames. We chose to rotate all frames to the same 
orientation before reducing them, so there is one extra 
interpolation involved in the processing of these frames, as 
the angles were such that a simple flipping of the images in 
x, y or both was not sufficient. 

/ band frames taken with ALFOSC have pronounced 
fringes - interference patterns due to the OH emissions 
coming from the upper ionosphere. If the fringes are stable 
throughout the night then it would be sufficient to obtain 
a fringe map by stacking non-aligned dithered I band 
frames and then subtract this from each individual image 
during reductions. This technique did not always remove 
the fringe pattern to our satisfaction, and residual "hills" 
and "valleys" could be clearly seen in the de-fringed images 
upon examination. Therefore we executed the pipeline in 
"optical pair-subtraction mode", using either dither or 
chopping mode pair-subtraction, as the data allowed. Since 
the fringes usually remain stable on such short time scales 
(~ 10 minutes from the start of one exposure to the next), 
this always produces clearly superior results compared to 
the average fringe map technique, and no fringe residuals 
can be seen on the de-fringed image even after smoothing 
it with a boxcar average (the smoothing is done solely for 



the purpose of inspection). 

EMMI There was a noticeable amplifier glow at varying 
positions along the edges of both CCDs in most / band 
frames of the 2005 data, which also seemed to change in 
time in position and intensity. If left uncorrected, this 
produces a significant non-flatness of the background in 
the stacked frame, which is impossible to remove with 
sky subtraction of a fitted plane. This issue is further 
enhanced by the fact that all four amplifiers were used 
during readout, each amplifier giving rise to its own varying 
amplifier glow. In such cases the pipeline was executed 
in "optical pair-subtraction mode", using either dither or 
chopping pair-subtraction, which successfully removed all 
amplifier glow. 

NOTCAM Approximately 8000 of the raw NOTCAM 
images (years 2002, 2004 and 2005) were lacking a world co- 
ordinate system (WCS) in the header. As the NIR pipeline 
relies on the WCS not only to align but also to identify 
2MASS stars for calibration, we had to find a way to 
introduce a WCS for every individual raw frame. Obtaining 
a full plate solution in any kind of semi-interactive fashion 
for this many images is not practical. Inste ad we used 
the p ublicly available astrometry.net software (|Lang et al.l 
|2010| ) to add the appropriate WCS to the headers. For 
most images the solver executed successfully using the 
standard indices built on the USNO-B catalog, but 
for a number of stubborn frames we successfully used 
2MASS indices instead. Since astrometry.net allows the 
user to search for solutions using already solved nearby 
fields, we were able to solve any remaining frames using 
local indices (e.g. an entire dither sequence of images 
of the same object can usually be solved as long as just 
one of those frames solves from USNO-B or 2MASS indices). 

The Wide-Field camera of NOTCAM suffers from an 
optical distortion which becomes severe towards the edges 
of the array, and can shift the projected location of sources 
by several arcseconds. We have applied a correction for 
this using the available distortion map on the NOTCAM 
website. This map is from the year 2005 and it worked well 
for our 2002, 2004, 2005 and 2006 NOTCAM data. We were 
unable to successfully remove all distortion for 2007 data 
with the same distortion map. For that data we chose to 
instead align the images using only the central 160 x 160 
arcsec (~ 700 x 700 pixels) around the galaxy. Since all 
targets obtained in 2007 have a fairly small apparent size, 
the lingering distortion along the edges is of no consequence 
for our analysis. 

SOFI There were no significant problems during the 
reductions of these data. All necessary calibration frames 
were of good quality, and a WCS was present in all frames. 

FORS The FORS exposures were taken for the purpose of 
re-calibrating existing observations, or to fill in for a num- 
ber of filters that were missing in our table of observations. 
There were no significant problems during the reductions. 
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Figure 2. A typical radial error composition (here for ES0185- 
13). The standard deviation of the mean inside a ring (o- sdom , 
open circles) dominates the error budget in the central regions of 
the galaxy. As the signal to noise ratio decreases with increasing 
radius the uncertainty in the sky (a s k y , filled squares) becomes the 
dominant error source. The composite error (solid line) is included 
in the errorbars of all surface brightness and color profiles. 



2.2 Photometric calibration 

All data w er e cali brated in the Veg a photometric s ystem . 
iHunt et alJ l| 19981 . Northern) and IPersson et all l|l99Sl . 
Southern) photometric standards were obtained for the 
NIR regime. However, we decided early on that the pipelines 
should calibrate the NIR data using secondary 2MASS 
standards on each raw frame. This way the photometric 
accuracy of the stacked NIR frames is nearly independent of 
the actual photometric conditions during the observations. 
For some targets the field of view contained too few visible 
stars, and calibration of each individual frame before 
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Figure 3. ALFOSC PSF obtained from stars in our frames. The 
PSF is the average of three stars, where one saturated star was 
used to obtain the far wings of the profile beyond ~ 12 arcsec. The 
individual star profiles were normalized to the total flux inside a 
5 arcsec radius before averaging. 



stacking was then impossible. In these cases a single zero 
point per galaxy per night was obtained and applied to 
the reduced stacked images since those always have more 
sources than the single exposure frames. 

In the optical the data were calibrated against images of 
Landolt standard stars obtained during each observing 
night. Occasionally we ran into problems when calibrating. 
Some nights were clear but non-photometric, others allowed 
for observations only in the beginning and the end of 
the night, resulting in an unacceptably long time interval 
between galaxy and standard star observations, and yet 
others had too few or no standard stars available because 
of saturation or bad pixels. Data from nights with bad 
atmospheric conditions were re-calibrated using short 1 or 
3 minute exposures, taken during photometric nights with 
FORS1 (VLT) or ALFOSC (NOT) in all cases where the 
zero point reliability was questionable. 

As we aim to reach very faint surface brightness levels, 
we felt it necessary to perform a sanity check on our 
calibration. We therefore compared the photometry of 
the stars in our frames to that of values taken from the 
Pickles stellar library in various color-color plots, both in 
the optical and NIR. This comparison showed no significant 
offset between the two, and we are therefore confident that 
our photometry is reliabkfl 

We estimated the zero point uncertainty for the optical data 
using the scatter of measured zero point values for each 
night. For the NIR we instead measured the mean offset 
in the magnitudes of stars on our frames when compared to 
2MASS photometry. This uncertainty is represented by a zp 

(§ED. 



1 These plots are available on demand. Email GM. 
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3 SURFACE BRIGHTNESS PROFILES 

BCGs often have irregular morphologies and it is therefore 
not always straightforward to decide on how best to 
obtain a radial distribution of the surface brightness. If the 
galaxy looks regular then ellipses will trace the isophotes 
quite well. If that is not the case then there are several 
options, developed specifically through the consideration 
of highly irregu lar starbursting galax i es, for how to ap- 
proach that (e.g. Papaderos et al.|[l996l ; ICairos et alj|2001al : 



iPapaderos et al.l I2002T ) . The usual assumption is that even 
if the starburst component is quite irregular, the underlying 
old stellar population is well-approximated by ellipses. 
We found, however, that if one goes deep enough, even 
that assumption is not always valid. We therefore decided 
to obtain the surface brightness profiles in two distinct, 
non-overlapping ways. 

First, we looked at surface brightness as a function of radius, 
i.e. we performed the integration in concentric elliptical 
rings, where each new data point is obtained by stepping in 
radius along the semi-major axis. The geometric parameters 
for the ellipse were obtained by adopting the shape of the 
ellipse best fitting the 25 th isophote as given by IRAF 
ELLIPSE. The parameters (ellipticity, position angle, and 
(x,y) ellipse center) were kept constant at each integration 
step. The position angle is measured counterclockwise from 
North to East. The best fitting ellipse was obtained for 
each target in a reference bandpass and then applied to 
the rest of the filters of that target. This way the physical 
area we integrate over at each radial step is identical in 
each filter. The reference bandpass is always taken to be 
the bluest available filter after the U band, which usually 
means the B band, or V if B is missing. We do not use the 
U band as reference on purpose since the atmosphere has 
poor transmission at those wavelengths and as a rule the U 
band is never deep enough. The B band on the other hand 
has long exposure times and a background which seems to 
be more stable than any other band at our disposal, which 
makes it a natural choice as a reference. Throughout this 
paper a step size of 1 arcsec was adopted for every bandpass 
and every target. This stepsize was chosen because 1 arcsec 
is on the order of the average seeing for all observations, 
which is ~ 1.1 arcsec. The B band ESO338-04 and HE2-10 
observations are an exception with ~ 1.7 arcsec seeing. 
Note that due to the different pixel scales between galaxies 
sometimes a 1 arcsec step would result in a non-integer 
number of pixels, which is rounded down to an integer 
and hence the actual step size will be slightly less than an 
arcsecond. More precisely, data with a pixel scale of 0.288 
arcsec have an actual step size of 0.86 arcsec. For each 
galaxy the images in all filters were resampled to a common 
(NIR) pixel scale. 

We further obtained isophotal radial profiles of the galaxies. 
These are interesting to compare to the standard elliptical 
profiles because they are essentially independent of the 
location of the star-forming knots and the morphology 
of the galaxy at any surface brightness level. There can 
be no "bumps" or "plateau" structures in an isophotal 
profile since at each step it samples the target in a manner 
which always produces a monotonously declining brightness 



distribution. In contrast to the elliptical integration profiles, 
an isophotal profile shows radius as a function of surface 
brightness. One always starts at the brightest pixel, and 
then steps to fainter isophote levels in magnitude bins. At 
each bin the number of pixels with surface brightness inside 
the current magnitude range is identified. This area may 
have any shape and does not necessarily need to be circular, 
elliptical or even connected. One then sets the area equal to 
an equivalent circular area and thus obtains an equivalent 
radius for each magnitude bin. At fainter isophotes this 
procedure breaks down as one starts to run out of available 
pixels and the pixels actually included in a magnitude bin 
start to become dominated by the positive outliers of the 
background noise. At this point it is common practice to 
instead approximate the fainter isophotes with ellipses and 
present a composite isophotal/elliptical profile. We have 
chosen not to do this here at any point along the isophotal 
profile, firstly because we already present standard profiles 
with elliptical integration and so an isophotal/elliptical 
composite integration profile will not give us any new 
information at large radii, and secondly because even at 
the faintest isophotal levels we have targets that are very 
poorly approximated by ellipses, e.g. UM1 60 and IIZw40 
in this sample and a number of galaxies in iMicheva et al.l 
|2012aT ). Instead, we chose to truncate the profiles at 
brighter isophotes where we are still confident that they are 
reliable. A constant magnitude bin of 0.5 mag was chosen 
for all galaxies and all filters. As for the elliptical profiles 
each isophotal area within the bins was defined using the 
reference filter (usually the B band) and then applied to 
the rest of the filters. This means that the physical area 
included in each bin is identical across all filters. 

Radial color profiles were obtained from both the elliptical 
integration and isophotal surface brightness profiles. The 
profiles are truncated when the composite error becomes 
too large, however that limit varies from color to color. 
Figure [T] presents the elliptical and isophotal surface 
brightness and color profiles for each galaxy. The profiles 
are reliable down to (Ik < 23 mag arcsec -2 in the NIR, 
and fiB < 28 mag arcsec -2 in the optical. Judging by the 
behavior of the corresponding profile and its errors, these 
limits seem to fit most galaxies except ESO400-43 (fj,v < 27 
mag arcsec -2 ), ES0462-IG 020 (p B < 26.5 mag arcsec -2 ), 
IIZw40 (hb < 25.5 mag arcsec -2 ), UM133 (fj, K < 22.5 
mag arcsec -2 ), and UM417 ((ik < 22 mag arcsec -2 ). On 
the other hand for several galaxies the profiles seem to 
reach much deeper than our imposed limit, with /j,k < 24 
mag arcsec -2 and fis J$ 29 mag arcsec -2 . We will further 
analyze the reliability of t he data in ij | 3.3| Ga lactic ex- 
tinction correction following ISchlegel et al.l ( 19981 ) has been 
applied to all filters, however, internal extinction has not 
been taken into account. Further, none of the profiles have 
been inclination-corrected because we found it impossible 
to define an unambiguous disk for too many of the galaxies 
in the sample. 

For every galaxy and every filter we examined the flatness 
of the residual background by smoothing the images with 
a 3 x 3 boxcar average. This was done solely for the 
purpose of looking for residual large scale structures, while 
unsmoothed images were used for the derivation of the 
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surface brightness profiles. 

Before obtaining elliptical and isophotal profiles all contami- 
nating foreground and background objects were masked out. 
The masking process uses SExtractor output as an indica- 
tion of each object's shape and size, and then increases each 
reported size by a factor of 2.5. While this works in most 
cases in fully automatic mode, we visually examined the 
sizes and shapes of masks for each individual frame, and 
adjusted them where necessary. 

3.1 Point spread function 

To test whether the point spread function (PSF) of the in- 
struments can affect the light distribution of a galaxy at faint 
levels, we obtained the PSF for ALFOSC, presented in Fig- 
ure [3] The radial light distribution of three stars, two faint 
and one very bright and saturated, were obtained from a 
random frame. Each profile was normalized by the flux sum 
inside a five arcsec radius, which scales the profiles down to 
the same magnitude levels. The average of the three profiles 
was taken where they overlap, which was out to a radius 
of ~ 12 arcsec. The saturated star was naturally not used 
for the central region of the PSF. Beyond ~ 12 arcsec the 
wings are based solely on the one saturated star. We then 
scaled the PSF up to the corresponding luminosity for a few 
random galaxies (observed with ALFOSC) to see if the in- 
strument PSF has any influence on the surface brightness 
profiles at large radii. Our test showed that it does not. 
The PSF drops in brightness much faster than our galaxies, 
therefore any light we detect at faint isophotal levels is not 
a PSF effect. 



3.2 Error estimation 

A very det ailed description of our error estimation is 
presented in iMicheva et al.l l|2010t) . however, here we will 
briefly outline how the errors are calculated. At each data 
point in the elliptical integration profiles we have accounted 
for 3 sources of uncertainty: a zp - the (constant) zero-point 
uncertainty, a s k y - the uncertainty in the flatness of the 
sky, and a a d om - the uncertainty in the mean intensity level 
of each elliptical ring, where "sdom" stands for standard 
deviation of the mean. As mentioned in § 12.21 the zero-point 
uncertainty is the mean deviation in the magnitude of all 
stars in our frames from their corresponding tabulated 
values in the Sloan Digital Sky Survey (SDSS) and 2MASS 
catalogs. If the galaxy is not found in the SDSS, then 
this uncertainty is instead the standard deviation of the 
photometric zero points obtained during the correspond- 
ing observing night. The sky uncertainty is obtained by 
calculating the mean intensity inside of square apertures, 
placed on "empty" sky regions in the frame. The standard 
deviation of these mean intensities then estimates the 
departure from the flat zero sky and is assumed to be 
constant in each pixel. The size of the apertures does 
matte r and we refer the interested reader to IMicheva et all 
(|2010h - Since we cannot measure where the sky level should 
truly lie, presenting a cr s ky which estimates the precision 
of the assumed sky level is meaningless and misleading. 
Therefore, our o s k y is designed to only be a measure of the 



flatness of the sky and also directly carries information on 
when the assumed sky would start influencing the signal 
from a target. However, this should not be interpreted as 
an estimate of the accuracy of the sky level, as it carries no 
information on whether or not the sky level in our images 
is true or not. The final source of uncertainty is in the 
accuracy with which we can measure the mean intensity 
inside each elliptical ring. This in itself is a composite error 
of the Poisson noise and the statistical fluctuations in flux 
across a ring and is given by the standard deviation of the 
mean (sdom) of each ring. It is a measure of the precision 
with which the mean intensity is known. This is not to be 
confused with the standard deviation inside each elliptical 
ring which is a measure of the scatter in flux - a physical 
property, not an error source. An example radial error 
composition for elliptical integration profiles is shown in 
Figure [2] The sky uncertainty clearly dominates the error 
budget at faint intensity levels. 

In an isophotal profile only the zero point and the sky back- 
ground uncertainties are accounted for. This is because by 
design here the standard deviation of the mean at each ra- 
dius will always be small and pretty much constant - pixels 
that fall inside the same magnitude bin all have very similar 
intensities so any small scatter they might have would not 
be a measure of anything useful. 

3.3 Systematic errors 

In IMicheva et all l|2010t ) we performed a test on the reliabil- 
ity of the reduction pipeline using synthetic disk galaxies. 
Specifically, we investigated at what surface brightness lev- 
els the systematic errors in the sky subtraction start dom- 
inating the observed flux. Since the systematic errors are 
generally unknown, this test allowed us to get a handle on 
their effects. We obtained a K band limit of \ik = 23 mag 
arcsec -2 and fiv = 28.5 mag arcsec -2 , beyond which the 
surface brightness profile would show a systematic deviation 
from the synthetic e xponential prof il e. We use the same re- 
duction pipeline as IMicheva et al.l (|2010h but these limits 
cannot be simply carried over to other data, as they are 
dependent on the depth of each observation, the individual 
accuracy of the natfielding, the typical behavior of the sky, 
etc. Regardless of the individual errorbars of each profile 
point if the data values appear to decrease/increase in an 
ordered fashion, then there is likely still useful signal mea- 
sured at the corresponding radius. If one is instead sampling 
the sky, then the data values should reflect random statis- 
tical fluctuations. Therefore one could naively expect that 
any ordered behavior would indicate the continued presence 
of a signal from the target galaxy. We choose to test this ex- 
plicitly for the B and K bands from the ALFOCS and NOT- 
CAM instru ments. Such a test a lready exists for the SOFI 
instrument (|Micheva et aljfeoiOj ). while the dither pattern 
on EMMI data left no space around the primary galaxy for 
a synthetic disk to be inserted and visible in each frame. 

In raw B band frames we inserted synthetic disks at six 
different positions (not simultaneously), after having "de- 
flattened" the disks by multiplying with the corresponding 
flatfield before insertion. The reduction and calibration of 
these frames then proceeded as usual. We obtained a surface 
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Figure 4. B band estimation of systematic errors. Upper 
panel: The surface brightness difference between a noiseless an- 
alytic exponential disk and the profile of the latter after a run 
through the reduction pipeline. The profiles of six exponential 
disks, placed at different positions in the raw frames, were ob- 
tained and averaged, and then subtracted from the analytic disk 
profile to obtain AfiB ■ Lower panel: Same as above for three disk 
positions, each with four different simulated flatfielding errors of 
up to 2% in amplitude. The errorbars are the standard deviation 
between the different profiles. 
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Figure 5. K band estimation of systematic errors. Upper 
panel: Same idea as in Figure [4] but for three exponential disks 
in the K band. Lower panel: Same as above for the three disk 
positions, but each with two different simulated flatfielding errors 
of up to 5% in amplitude. The errorbar budget is the same as in 
Figure [4] 



brightness profile for each disk position and compared it 
to the profile of the analytical exponential disk that was 
used. In all cases we recover the original properties of the 
disk - scale length, central surface brightness, Holmberg 
radius, total luminosity down to the disk truncation - to 
within 0.05 — 0.1 units. In the top panel of Figure [4] we 
show the magnitude difference between the original analytic 
exponential disk and the recovered surface brightness 
profile, where the latter is the average of the six different 
cases. The plot shows that there is no difference between 
the two down to hb ~ 26 mag arcsec -2 , and a systematic 
brightening of the averaged recovered profile beyond this 
isophote, which is on average smaller than A/i_b ~ 0.1 mag 
down to lib ~ 28.5 mag arcsec -2 . Judging from the plot, 
the profile is probably still reliable even down to \ib ~ 29 
mag arcsec -2 . To test the effect of enhanced flatfielding 
errors, for three of the positions we also simulated four 
cases of flatfielding errors of varying amplitude but with a 
maximum deviation of 2%, and then proceeded as described 



Table 3. Surface photometry parameters. The position angle 
(PA°) in degrees and the cllipticity (e) were obtained from the B 
band image and applied to the remaining filters. PA° is measured 
counterclockwise from North to East. The radius where the mean 
surface brightness is fiB ~ 26.5 mag arcsec -2 is the Holmberg 
radius, Rh, and is given in arcseconds and kpc. The absolute B 
magnitude, Mb, is calculated from th e area inside Rh an d has 
been corrected for Galactic extinction l lSchleeel et al.ll 199811 . The 
Mb errors are identical to the errors of the apparent B magnitude 
in Tabled 



Galaxy 


PA° 


e 


Rlj 


R k * c 

ri 


M B 


ES0185-13 


46 


0.18 


22.5 


8.3 


-19.48 


ES0249-31 


45 


0.39 


129.6 


6.8 


-17.99 


ESO338-04 


76 


0.45 


43.8 


7.9 


-19.31 


ESO400-43 t 


-15 


0.22 


23.9 


9.2 


-20.15 


FSO421-02 


—89 


0.26 


56.2 


3.4 


— 16.53 


ESO462-20 


-86 


0.40 


27.6 


10.6 


-20.03 


HE2-10 


28 


0.15 


96.8 


7.4 


-19.07 


HL293B 


-2 


0.12 


11.2 


0.9 


-13.71 


IIZW40 


-71 


0.19 


109.7 


6.3 


-18.68 


MK600 


-51 


0.48 


44.9 


2.4 


-15.12 


MK900 


41 


0.27 


49.6 


2.7 


-15.99 


MK930 


41 


0.19 


19.7 


6.8 


-19.23 


MK996 


20 


0.29 


27.6 


2.4 


-16.25 


SBS0335-052E 


-61 


0.11 


11.2 


2.9 


-16.89 


SBS0335-052W 


-50 


0.22 


8.6 


2.3 


-14.53 


TOL0341-407 


-86 


0.40 


19.0 


5.6 


-17.50 


TOL1457-262 I 


-21 


0.54 


32.8 


11.6 


-20.02 


TOL1457-262 II 


85 


0.74 


56.2 


19.8 


-19.54 


UM133 


16 


0.40 


37.4 


3.3 


-15.81 


UM160 


SO 


0.14 


48.4 


6.6 


-17.76 


UM238 


-78 


0.51 


22.5 


5.9 


-16.76 


UM417 


-17 


0.40 


14.0 


2.3 


-14.50 


UM448 


11 


0.24 


43.8 


17.5 


-20.42 


UM619 


45 


0.46 


26.8 


8.9 


-18.40 



f - V instead of B filter 

above to obtain the magnitude difference as a function 
of surface brightness level. This is presented in the lower 
panel of Figure [4] We reach lib ~ 28.0 mag arcsec -2 with 
the same excess of 0.1 on average, and /is ~ 29.0 mag 
arcsec - with an average of 0.2 mag offset. In light of these 
results, we are confident that we can trust the errorbars 
in our surface brightness profiles for the sample data, and 
we set our reliability limit in the B band to a conservative 
(ib ~ 28.0 mag arcsec -2 . 

For the K band in a similar fashion we find negligible A/hk 
down to (j,k ~ 22.5 mag arcsec -2 (Figure [5}, and an offset 
of max 0.2 mag down to (xk ~ 23.3 mag arcsec -2 . Therefore 
we set the reliability limit in the K band to lik ~ 23.0 mag 

arcsec -2 . 



4 CONTOUR PLOTS 

The isophotal contour plots for each galaxy were obtained 
from the deepest best image (the B band for the majority 
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Table 4. Integrated surface photometry for the sample. The integration is carried out down to the Holmberg radius R H , which 
is de fined from jig for ea ch target and then applied to the remaining filters. All values have been corrected for Galactic extinc- 
tion jSchleeel et al.lll99gfi . 



Galaxy 


B 


U — B 


B — V 


V — R 


V — I 


V - K 


H — K 


ES0185-13 


14.94 ± 0.05 


-0.17 ± 0.16 


0.26 ± 0.05 


0.19 ± 0.03 


0.49 ± 0.03 


1.89 ± 0.06 


0.19 ± 0.06 


ES0249-31 


12.18 ± 0.04 




0.50 ± 0.05 


0.31 ± 0.03 


0.66 ± 0.03 


2.12 ± 0.04 


0.15 ± 0.04 


ESO338-04 


13.55 ± 0.04 




0.34 ± 0.05 


0.07 ±0.05 


-0.22 ± 0.06 


1.04 ±0.06 


0.11 ±0.05 


ESO400-43 


14.34 ± 0.19* 








0.38 ± 0.20 


1.83 ± 0.20 


0.25 ± 0.05 


ESO421-02 


13.94 ± 0.04 


-0.52 ± 0.15 


-0.15 ± 0.05 


0.63 ± 0.03 


0.97 ±0.03 


2.66 ±0.07 


0.19 ±0.07 


ESO462-20 


14.47 ± 0.04 


-0.16 ± 0.16 


0.29 ± 0.05 






2.03 ± 0.19 


0.44 ± 0.20 


HE2-10 


11.93 ± 0.04 


0.11 ± 0.14 


0.48 ±0.05 






2.78 ±0.06 


0.45 ±0.08 


HL293B 


17.35 ± 0.05 


-0.50 ± 0.16 


0.48 ± 0.07 


0.11 ± 0.07 


0.22 ± 0.09 


1.60 ±0.08 


0.19 ± 0.15 


IIZW40 


11.68 ± 0.04 


-0.94 ± 0.18 


0.22 ± 0.06 




-0.36 ± 0.09 


0.97 ± 0.07 


1.25 ± 0.09 


MK600 


15.07 ± 0.06 


-0.31 ± 0.16 


0.37 ± 0.07 


0.20 ± 0.07 


0.41 ± 0.10 


1.38 ±0.05 


-0.10 ± 0.11 


MK900 


14.25 ± 0.04 


-0.20 ±0.15 


0.58 ± 0.06 


0.44 ± 0.07 


0.81 ±0.09 


2.34 ±0.10 


0.20 ±0.10 


MK930 


15.04 ±0.05 


-0.77 ±0.16 


0.41 ± 0.06 


0.07 ±0.07 


0.17 ±0.09 


1.86 ±0.09 


0.42 ±0.09 


MK996 


15.05 ± 0.02 


-0.30 ±0.15 


0.57 ±0.05 


0.42 ± 0.06 


0.78 ± 0.09 


2.47 ±0.05 


0.28 ±0.09 


SBS0335-052 E 


16.77 ±0.05 








-0.02 ±0.10* 


1.19 ± 0.05* 


0.72 ±0.07 


SBS0335-052 W 


19.14 ±0.05 








0.29 ±0.09* 


1.62 ± 0.44* 


0.22 ±0.84 


TOL0341-407 


16.41 ± 0.04 


-0.51 ± 0.18 


0.66 ± 0.05 




0.20 ±0.50 


1.63 ±0.06 


0.17 ±0.07 


TOL1457-262 I 


14.29 ± 0.05 




0.28 ± 0.06 






1.61 ±0.07 




TOL1457-262 II 


14.78 ± 0.20 




0.19 ± 0.29 






2.25 ±0.21 




UM133 


15.50 ±0.05 


-0.43 ± 0.16 


0.34 ±0.06 


0.30 ±0.07 


0.58 ±0.10 


2.14 ±0.09 


0.34 ±0.10 


UM160 


14.47 ±0.02 


-0.95 ± 0.06 


0.42 ± 0.04 


0.16 ±0.05 


0.47 ±0.05 


1.83 ±0.16 


0.13 ±0.16 


UM238 


16.91 ± 0.05 


-0.33 ± 0.16 


0.40 ±0.07 


0.30 ±0.07 


0.62 ±0.10 


2.12 ±0.16 


0.57 ±0.19 


UM417 


18.15 ± 0.02 


-0.44 ±0.04 


0.58 ± 0.04 


0.13 ±0.04 


0.22 ±0.06 


1.82 ±0.45 


0.98 ±0.48 


UM448 


14.16 ± 0.02 




0.49 ± 0.04 


0.26 ±0.06 


0.64 ±0.05 


2.37 ±0.12 




UM619 


15.78 ± 0.03 


-0.10 ±0.09 


0.38 ± 0.05 


0.36 ±0.06 


0.69 ±0.06 


2.59 ±0.20 


0.81 ±0.24 



f - V instead of B filter 
\ - B instead of V filter 



of the cases). The python astLib package was used to re- 
sample the images to tangential projection so that North 
is up and East is to the left, and to obtain the axis labels 
in terms of RA and Dec, but the plotting of the contours 
was done with the built-in pylab function 'contour' in bins 
of 0.5 magnitudes. Beyond a certain surface brightness the 
isophotes become noisy as the magnitude bins start sampling 
the sky background. Therefore beyond this level, which is 
individually decided upon for each galaxy and ranges from 
(is ~ 22.0 to /is ~ 27.5 mag arcsec -2 , the images were 
smoothed with a boxcar median filter of increasing size (e.g. 
5x5, 15x15, 21x21, etc.) for each consecutive magni- 
tude bin (e.g. /j, B = 24 - 24.5, 24.5 - 25, 25 - 25.5 mag 
arcsec -2 , etc.). The upper right panel of Figure [T] shows B 
band contour plots for each galaxy. 



5 RGB IMAGES 

RGB images give an invaluable overview of the morphology 
of a galaxy, and we provide such images for each galaxy with 
the three channels (Blue, Green, Red) assigned to (U, B, I) 
in the majority of the cases. We tested various available 
tools for scaling and combining the RGB channels, includ- 
ing GIMP, FITS Liberator with GIMP, DS9, Aladvn, and 
the python astronomical plotting package APLpy. In most 
of these cases we find that it is very hard to select channel 
scaling factors that would universally give good results for 
all galaxies. In addition some of these algorithms regularly 
saturate the color of a channel in high-intensity areas of the 



image, thus obscuring all features. The best re sults were ob- 
tained with our own implementation of the iLupton et al.l 
(2004) algorithm. Only APLpy gave results of compara- 
ble quality but required channel adjustments for individual 
galaxies. In order to facilitate direct comparison one must 
avoid having di fferent channe l scalin g factors between galax- 
ies, which the ILupton et alj (I2004T ) algorithm allows. The 
benefit of this algorithm is that the intensity of a channel 
may saturate in certain parts of the galaxy, but the color 
never does, thus preserving all interesting morphological fea- 
tures that some of the more irregular galaxies in our sample 
have. The RGB channels were corrected for Galactic ex- 
tinction and converted to the AB magnitude system before 
scaling and combining. The lower right panel of Figure [T] 
shows the RGB image of each galaxy. To adjust the readers 
eyes to our color-scheme, in Figure [5] we show an example 
of SDSS multicolor images and our own. 



6 INTEGRATED SURFACE PHOTOMETRY 

TableUsummarizes the results of the integrated surface pho- 
tometry for all targets. The total luminosity of a galaxy was 
obtained by summing the flux inside of circular/elliptical 
apertures up to the Holmberg radius Rh- For all targets ex- 
cept for ESO400-43 the B band was used to define the Holm- 
berg radius and the shape and placement of the aperture, 
which was then applied to all remaining filters. The abso- 
lute magnitude was obtained with the cosmology-corrected 
(Hq = 73 km~ Y s~ x Mpc~ x , Sl ma ttcr = 0.27) Heliocentric dis- 
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Figure 6. UM160 (top) and UM448 (bottom) in SDSS (left) 
and our (right) RGB color schemes. 



tances available at the NASA/IPAC Extragalactic Database 
(NED) for each target. We are missing the B filter for 
ESO400-43, so the Holmberg radius and the total luminos- 
ity were instead obtained from the V band. Similarly, for 
SBS0335-052E&W we present B - I and B - K photome- 
try due to lack of V band. The presented magnitudes have 
not been inclination-corrected because the perturbed mor- 
phology of some targets does not allow for an unambiguous 
identification of a disk structure. 

6.1 Error estimation 

The errors for the total luminosities were estimated by vary- 
ing the geometric parameters (ellipticity and position angle) 
of the aperture around the best-fitting values used in the 
surface brightness profile integration, listed in Table 2] The 
position angle was varied by ±5°, the ellipticity by ±0.1. 
The standard deviation of the luminosity values obtained in 
such a fashion was added in quadrature to the photometric 
uncertainty for each filter and target. No other errors were 
included since in this case the sky uncertainty cannot signif- 
icantly influence the flux sum of the entire galaxy, and the 
flux fluctuations inside the aperture reflect a physical prop- 
erty of the galaxy (radially changing morphology), rather 
than being a measure of uncertainty. The errors in the total 
apparent and total absolute magnitudes are the same since 
no distance errors were included in the error budget. 



7 TOTAL COLORS 

In order to interpret the observed total colors of the 
galaxies, we must first separate the stellar populations. In 
a BCG where the old population displays regular isophotes 



one can easily model it with a disk and subtract it in order 
to obtain the flux from the young population, regardless of 
how concentrated or extended the latter may be. Most of 
the galaxies in our sample, however, do not have regular 
isophotes even in the outskirts, and hence it would be 
difficult to estimate the reliability of such a subtraction. We 
further wish to subject all galaxies to the same treatment, so 
that any differences in results will unambiguously indicate 
differences in galaxy properties and morphology, and will 
not be simply the result of having individual definitions 
for each galaxy. However, upon examining the individual 
surface brightness profiles for all filters and galaxies, it 
became apparent that using one single definition for what 
comprises the outskirt regions of a galaxy, or how extended 
the starburst is, leads to less than informative results for 
some particularly troublesome targets. Furthermore, the 
sky brightness in the NIR images is higher than in the 
optical, and this limits the radial range to which we can 
trace the profiles in all filters simultaneously. For example, 
a change in population is clearly visible for ES0185-13 
at fiB ~ 26.5 mag arcsec -2 and radius r ~ 20 arcsec, 
but at that radius either of the NIR profiles is completely 
dominated by sky fluctuations. For other targets, e.g. 
Toll457-262/&77, the NIR profiles are still reliably trace- 
able at radii corresponding to /is ~ 26 - 28 mag arcsec -2 . 
Despite the deep NIR images, which have systematically 
longer exposure times than in the optical by a factor of 
three on average, it is obviously impossible to both study 
the old extended population beyond /is ~ 26 mag arcsec -2 
and simultaneously only use radial ranges with reliable NIR 
data. This is due to the fact that ground-based observations 
in the NIR are not photon-limited and increasing the 
exposure time cannot push down the detected isophotes to 
arbitrarily faint levels. Therefore we divide our investigation 
in separate radial ranges: the center region is defined as the 
ellipse whose semi-major axis is the radius corresponding 
to /is ~ 24 mag arcsec -2 ; the intermediate (young and 
old population mix) region as the elliptical ring between 
24 < Ms < 26 mag arcsec -2 ; and the galaxy outskirts 
region as between 26 < /is < 28 mag arcsec -2 . The ellipse 
parameters used during this integration are the same as in 
Table 3] For color-color plots involving purely optical colors 
(U — B, B — V, V — R, V — I) we can safely use data from 
the outermost region. For plots involving optical-NIR color 
combinations (V — K, B — K) we restrict ourselves to using 
the intermediate region, even though it is often unclear 
whether the burst contribution in that region is negligible. 



The total colors for these radial ranges are presented in 
Table [5] Each value was calculated by summing up all pixels 
inside the region, after foreground stars and background 
galaxies have been masked out. As already mentioned 
above, each region was defined from a reference filter, 
usually the B band, and then applied to the rest of the 
filters. Additionally, a union of the source masks in all filters 
was created and applied to the region in each filter before 
summation. Note that this means that the summation is 
carried over the same physical region of the galaxy in each 
filter and involves the exact same number of pixels. While 
it pains us to throw away perfectly good pixels this way, 
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Table 5. Total colors for radial ranges corresponding to fj,g < 24 , 24 < , lib < 26, an d 26 < fiB ;$ 28 mag 
arcsec - 2 . All values have been corrected for Galactic extinction l lSchlegel et al.|[l998l) . The errors include 
°"sfcyi a sdom and &zp added in quadrature. 



Gl alaxy 






U - B 


B - V 


V — R 


V - I 




V - K 


H - K 


ES0185- 


13 


24 
26 


-24 
-26 
-28 


-0.16 ± 0.16 
-0.39 ± 0.19 
-0.76 ± 0.51 


0.26 ± 0.06 
0.31 ± 0.10 
0.29 ± 0.47 


0.15 ± 0.05 
0.52 ± 0.10 
0.72 ± 0.50 


0.45 ± 0.05 
0.81 ± 0.12 
0.40 ± 0.88 


1 
1 


.87 ± 0.07 
.93 ± 0.67 


0.21 ± 0.07 
-0.13 ± 0.94 


ES0249- 


31 


2-1- 
26 


-24 
-26 
-28 




0.47 ± 0.05 
0.59 ± 0.05 
0.61 ± 0.13 


0.34 ± 0.03 
0.30 ± 0.03 
0.08 ± 0.16 


0.65 ± 0.03 
0.74 ± 0.04 
0.70 ± 0.19 


2 
1 


.22 ± 0.05 
.64 ± 0.31 


0.20 ± 0.04 
-0.25 ± 0.35 


ES0338- 


04 


24- 
26 


-24 
-26 
-28 




0.30 ± 0.09 
0.27 ± 0.12 
0.10 ± 0.59 


0.04 ± 0.06 
0.36 ± 0.09 
0.22 ± 0.66 


-0.27 ± 0.07 
0.66 ± 0.11 
1.06 ± 0.58 


0.99 ± 0.06 
1.62 ± 0.41 


0.09 ± 0.07 
-0.12 ± 0.78 


ESO400- 


43 


*-24* 
24-26* 
26-28* 








0.37 ± 0.20 
0.58 ± 0.22 
0.75 ± 0.65 


1 
1 


.85 ± 0.20 
.51 ± 0.50 


0.28 ± 0.07 
-0.31 ± 0.59 


ES0421- 


02 


24 
26 


-24 
-26 
-28 


-0.54 ± 0.15 
-0.35 ± 0.17 
-0.79 ± 0.36 


-0.20 ± 0.06 
0.08 ± 0.11 
0.41 ± 0.43 


0.63 ± 0.04 
0.64 ± 0.11 
0.36 ± 0.47 


0.94 ± 0.05 
1.08 ± 0.11 
1.07 ± 0.45 


2 
2 
2 


.67 ± 0.08 
.72 ± 0.16 
.69 ± 0.66 


0.18 ± 0.08 
0.25 ± 0.20 
0.63 ± 1.04 


ES0462- 


20 


24 
26 


-24 
-26 
-28 


-0.11 ± 0.16 
-0.52 it 0.26 
-1.22 ± 1.75 


0.26 ± 0.07 
0.63 ± 0.21 
1.58 ± 1.72 






2 
1 


.06 ± 0.20 
.41 ± 0.48 


0.48 ± 0.20 
-0.08 ± 1.01 


HE2-10 




24 
26 


-24 
-26 
-28 


0.05 ± 0.14 
0.26 ± 0.24 
0.27 ± 0.88 


0.41 ± 0.06 
0.68 ± 0.11 
0.62 ± 0.45 






2 
2 
2 


.85 ± 0.07 
.65 ± 0.30 
.38 ± 1.67 


0.46 ± 0.09 
0.50 ± 0.41 
0.45 ± 2.30 


HL293B 




24 
26 


-24 
-26 
-28 


-0.52 ± 0.17 
-0.46 ± 0.17 
-0.39 ± 0.37 


0.47 ± 0.09 
0.51 ± 0.08 
0.69 ± 0.20 


0.07 ± 0.10 
0.33 ± 0.09 
0.28 ± 0.24 


0.14 ± 0.11 
0.57 ± 0.12 
0.60 ± 0.28 


1 

2 


.43 ± 0.15 
.25 ± 0.90 


0.24 ± 0.19 
0.09 ±1.16 


IIZW40 




24 
26 


-24 
-26 
-28 


-0.82 ± 0.23 


0.18 ± 0.07 
0.34 ± 0.74 
0.07 ± 2.61 




-0.27 ± 0.11 


0.89 ± 0.09 
1.23 ± 1.26 


0.34 ± 0.31 


MK600 




24- 
26 


-24 
-26 
-28 


-0.36 ± 0.16 
-0.12 ± 0.21 
-0.14 ± 0.83 


0.37 ± 0.07 
0.39 ± 0.09 
0.34 ± 0.42 


0.16 ± 0.07 
0.29 ± 0.09 
-0.26 ± 0.50 


0.37 ± 0.09 
0.59 ± 0.14 
-0.84 ± 2.16 


1 
1 


.60 ± 0.09 
.14 ± 1.99 


0.20 ± 0.12 
-0.56 ± 2.15 


MK900 




24 
26 


-24 
-26 
-28 


-0.24 ± 0.16 
0.02 ± 0.17 
-0.02 ± 0.48 


0.55 ± 0.07 
0.70 ± 0.07 
0.66 ± 0.20 


0.43 ± 0.07 
0.51 ± 0.08 
0.49 ± 0.23 


0.76 ± 0.09 
0.92 ± 0.10 
1.01 ± 0.26 


2 
2 
2 


.26 ± 0.11 
.54 ± 0.27 
.67 ± 1.45 


0.18 ± 0.10 
0.22 ± 0.34 
0.43 ± 2.05 


MK930 




24 
26 


-24 
-26 
-28 


-0.79 ± 0.16 
-0.68 ± 0.18 
-0.31 ± 0.55 


0.43 ± 0.08 
0.28 ± 0.12 
-0.03 ± 0.67 


0.05 ± 0.08 
0.39 ± 0.13 
0.15 ± 0.92 


0.14 ± 0.11 
0.58 ± 0.16 
-0.03 ± 1.41 


1 
2 


.82 ± 0.10 
.44 ± 0.72 


0.48 ± 0.10 
-0.00 ± 0.85 


MK996 




24 
26- 


-24 
-26 
-28 


-0.32 ± 0.16 
-0.09 ± 0.16 
-0.28 ± 0.32 


0.56 ± 0.06 
0.65 ± 0.05 
0.59 ± 0.16 


0.42 ± 0.07 
0.47 ± 0.07 
0.47 ± 0.16 


0.77 ± 0.10 
0.91 ± 0.10 
0.88 ± 0.23 


2 
2 
2 


.48 ± 0.07 
.44 ± 0.15 
.13 ± 1.33 


0.32 ± 0.11 
0.15 ± 0.52 
-1.02 ± 2.01 


SBS0335 


-052E 


24 


-26 








-0.07 ± 0.15* 
0.30 ± 0.18* 


1.10 ± 0.14* 
1.99 ± 0.86* 


0.54 ± 0.17 






26 


-28 








0.65 ± 0.74* 








SBS0335 


-052W 


24 
26 


-24 
-26 
-28 








0.18 ± 0.13* 
0.30 ± 0.20* 
0.60 : 0.64* 


1.32 ± 0.22* 
1.97 ± 0.77* 
2.86 ± 1.81* 


0.53 ± 0.42 
0.15 ± 1.25 
-0.03 ± 2.63 


TOL0341-407 


2-1- 
26 


-24 
-26 
-28 


-0.51 ± 0.16 
-0.57 ± 0.17 
-0.53 ± 0.55 


0.67 ± 0.06 
0.67 ± 0.07 
0.72 ± 0.30 




0.14 ± 0.50 
0.55 ± 0.50 
0.89 ± 0.60 


1 
2 


.63 ± 0.09 
.01 ± 0.42 


0.23 ± 0.09 
0.25 ± 0.46 


TOL1457-262/ 


24 
26 


-24 
-26 
-28 




0.28 ± 0.07 
0.27 ± 0.10 
0.24 ± 0.48 






1 
2 


.61 ± 0.07 
.77 ± 0.22 




TOL1457-262/7 


24 
26 


-24 
-26 
-28 




0.19 ± 0.06 
0.23 ± 0.09 
0.35 ± 0.37 






2 
2 


.22 ± 0.06 
.56 ± 0.19 




UM133 




24 
26 


-24 
-26 
-28 


-0.48 ± 0.16 
-0.25 ± 0.16 
-0.05 ± 0.39 


0.30 ± 0.06 
0.40 ± 0.07 
0.39 ± 0.18 


0.24 ± 0.06 
0.36 ± 0.07 
0.37 ± 0.23 


0.48 ± 0.09 
0.75 ± 0.10 
0.97 ± 0.27 


1 
2 


.95 ± 0.23 
.33 ± 0.95 


0.30 ± 0.28 
0.36 ± 1.16 


UM160 




24 
26- 


-24 
-26 
-28 


-0.98 ± 0.06 
-0.89 ± 0.09 
-0.76 ± 0.59 


0.45 ± 0.04 
0.40 ± 0.06 
0.47 ± 0.34 


0.21 ± 0.05 
0.19 ± 0.09 
0.18 ± 0.54 


0.55 ± 0.05 
0.44 ± 0.11 
-0.01 ± 1.06 


1 
1 
2 


.92 ± 0.13 
.99 ± 0.56 
.78 ± 1.93 


0.15 ± 0.13 
0.08 ± 0.61 
0.08 ± 2.11 


UM238 




24- 
26 


-24 
-26 
-28 


-0.32 ± 0.16 
-0.30 ± 0.17 
-0.16 ± 0.44 


0.40 ± 0.07 
0.40 ± 0.07 
0.40 ± 0.22 


0.32 ± 0.07 
0.30 ± 0.08 
0.04 ± 0.34 


0.64 ± 0.10 
0.54 ± 0.14 
-0.52 ± 1.40 


2 
1 
2 


.13 ± 0.10 
.98 ± 0.57 
.82 ± 1.38 


0.55 ± 0.16 
0.46 ± 0.99 


UM417 




24 
26 


-24 
-26 
-28 


-0.58 ± 0.09 
-0.28 ± 0.08 
0.05 ± 0.54 


0.62 ± 0.09 
0.42 ± 0.05 
0.66 ± 0.21 


-0.08 ± 0.09 
0.29 ± 0.06 
0.23 ± 0.21 


-0.06 ± 0.09 
0.59 ± 0.08 
0.47 ± 0.34 


1 
2 


.48 ± 0.31 
.20 ± 0.75 


0.92 ± 0.66 
0.91 ± 1.67 


UM44S 




24- 
26 


-24 
-26 
-28 




0.47 ± 0.06 
0.58 ± 0.05 
0.74 ± 0.09 


0.24 ± 0.07 
0.41 ± 0.07 
0.33 ± 0.16 


0.60 ± 0.07 
0.93 ± 0.07 
0.96 ± 0.15 


2 
2 


.43 ± 0.14 
.09 ± 1.57 




UM619 




24- 
26 


-24 
-26 
-28 


-0.10 ± 0.10 
0.05 ± 0.18 
-0.13 ± 0.64 


0.36 ± 0.06 
0.44 ± 0.07 
0.59 ± 0.18 


0.34 ± 0.05 
0.59 ± 0.07 
0.82 ± 0.19 


0.68 ± 0.06 
0.92 ± 0.08 
1.17 ± 0.22 


2 
2 
3 


.31 ± 0.10 
.62 ± 0.64 
.05 ± 1.89 


0.52 ± 0.18 
0.89 ±1.13 





t - V instead of B filter 
X - B instead of V filter 
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possibly with high signal-to-noise, not doing so could bias 
the results and unnecessarily complicate the interpretation. 

At faint surface brightness levels the flux in some pixels will 
start randomly falling below a value of 0, which is an effect 
of the sky subtraction. To avoid a systematic reddening of 
colors, all pixels are summed, positive and negative, and no 
clipping of any kind is applied. Since the mean sky level is 
~ in each frame, the negative sky pixels cancel out the 
positive ones and we recover total colors consistent with the 
trends observed in the radial color profiles for a given magni- 
tude range. There is of course always a sky residual present 
in the frames, the total sum of the sky is never exactly re- 
gardless of what sky subtraction technique is used, but the 
influence of that sky residual is estimated and included in 
the total errors in the form of <7 s ky 

7.1 Error estimation 

The errors in the total colors were obtained in a fashion 
similar to £j |3.2l for each filter and then added in quadrature, 
<j(X — Y) = J a x + <7y , where ax and ay each include all 
three sources of uncertainty - <j s k y , &sdom and a zp - Note that 
only Osdom varies with each region (center, intermediate, and 
outskirts), while the other two are constant over the entire 
frame in each filter. 



8 HOST STRUCTURAL PARAMETERS 

The scale length h B and the central surface brightness n B 
are obtained for each galaxy by a weighted least squares fit 
of an exponential disk to the reference profile in the two 
radial ranges ARi and AR2 defined by (ib = 24-26 and 
(is = 26-28 mag arcsec -2 respectively. When applied to 
the NIR filters, the radial ranges ARi and especially AR2 
usually encompass a very noisy part of the profile. Even a 
weighted fit to such data points will give exponential disk pa- 
rameters with arbitrarily large uncertainties. In other words, 
we often cannot measure how the light distribution at NIR 
wavelengths behaves at such large radial distances from the 
center. Therefore fitting individual scale lengths for each fil- 
ter at the same AR\ and AR2 will give essentially random h r 
and (io- Since we do not know the true profile shapes at these 
radii, we make the assumption that the h B scale length, 
which is known with a fairly high accuracy as the B band 
profile is usually the deepest, is valid for all wavelengths and 
only fit the central surface brightness / jq in each filter. This 
assum ption is not unreasonable, e.g. iFerreiro fc Pastorizal 
l|2004h find that the disk scale length of minor mergers does 
not change in BVI. The same appears to be true for nor- 
mal ellipticals, as seen from the radial color profiles of the 
red and dead elliptical galaxy KIG 732 in Appendix [X] In 
hybrid optical-NIR colors strong gradients can still exist as 
demonstrated by the V — K radial color profile for KIG 732. 
This means that our burst estimate at NIR wavelengths is 
systematically overestimated. Further, in the NIR regime 
burst estimation often fails due to the complete breakdown 
of the NIR profile inside ARi or AR2, but it gives reason- 
able estimates of fio for profiles which still exhibit a statis- 
tically significant ordered signal, only one with very large 
errors. We will refrain from overanalysing the NIR burst 



colors, and will predominantly use optical burst estimates. 
Host structural parameters for all galaxies are presented in 
Table H 



8.1 Error estimation 

The h B errors are the propagated errors of the fitted slope, 
scaled to units of the profile errors at those radii. The /i B 
errors include the error of the fit and the zero point errors 
<j zp . The fio errors for all other filters also account for the 
unsealed h B error since the uncertainty in the assumed slope 
directly influences the extrapolated value. 



9 ESTIMATION OF THE BURST 
COMPONENT 

The central region of the profiles is usually dominated by 
the strong signal from the starburst population but it also 
contains a non-negligible contribution from the underlying 
old population. Notable exceptions to this are irregular mor- 
phology targets (e.g. Tol0341-407, UM160, UM238, etc.) for 
which the center of light (r = arcsec) is not at all located 
at the burst as there are several knots of star formation 
spread around the galaxy. If one can disentangle the contri- 
butions from the two populations one can obtain the colors 
of the burst. This is a rather complicated problem without 
a universal solution and lies beyond the scope of this paper. 
However, a simplified estimate of the light of the burst in 
each filter can be obtained by a ID subtraction of the fitted 
exponential disk from the profile data. For each filter, we use 
the B band h r and the fitted for the individual filters to 
interpolate along the exponential disk line, thus obtaining 
an estimate of the contribution of the old population to the 
total flux at that radius. The excess in luminosity above the 
exponential disk is an estimate of the burst luminosity. Rel- 
ative burst contribution to the B band and apparent burst 
luminosity obtained in such a fashion for all filters and galax- 
ies are presented in Table [7] We note that the morphology 
of some galaxies clearly contradicts the assumption of an 
underlying regular disk, even though an exponential profile 
can still be fitted to the light distribution at large radii. One 
should keep this in mind and refrain from overinterpreting 
the data. 



9.1 Error estimation 

The errors of the burst component parameters are estimated 
from the fitting errors of the exponential disk, scaled to units 
of the profile errors, and combined with the zero point un- 
certainty. Thus they may be underestimated since there is 
no measure of the uncertainty in the exact flux level. There 
is also no sky uncertainty contribution to the total error 
budget of the burst, but as Figure [2] shows, at such bright 
isophotal levels is it always the uncertainty in the flux level 
that would dominate. All errors are assumed to be indepen- 
dent and added in quadrature. 
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10 ASYMMETRY AND CONCENTRATION 

Asymmetry, dumpiness and concentration have been argued 
to contain information on the past and present evolution 
hist ory of a galaxy even when measured in only one fil- 
ter jConselicc 2003), which makes them uniquely economical 
parameters to obtain. The asymmetry is thought to corre- 
late with galaxy colors in a way which can separate re cent 
major mergers from non-mergers (|Conselice et al.l|2000h . To 
investigate this, in Table [5] we compute the asymmetries for 
all targets using the method 

El*>-J#| 



2E|/«I 



(1) 



because it correlates well with the colors, and hence is a 
good indicator of star formation CCons elice et al.l [2000). 
Io is the original (processed) image, and 1$ is rotated 180 
degrees. Decreasing signal-to-noise (S/N) in the fainter 
isophotal levels can greatly affect the minimum asymme- 
try, hence normally one would apply a noise correction. 
All galaxies in the sample are either starbursting or 
star-forming at moderately high SFRs, however, and the 
composite asymmetry due to both star formation (flocculent 
asymmetry) and merger (dynamical asymmetry) events is 
likely to be characterized by the presence of the bright star 
forming regions and not dominated by any possible faint 
extended tidal structures. Hence, instead of applying a 
noise correction to reach fainter isophotal levels, we choose 
to compute the asymmetry inside the Petrosian radius 
r[rj = 0.2], which typically gives asymme try errors lower 
than 0.02 rms for the galaxies in the sample (|Conselice et al.l 
2000). 77 denotes the inverted 77, defined as the ratio between 
the local surface brightness at some radius and the average 
surface brightness inside that radius. The location of the 
starburst regions in the individual targets is varied, and 
they can be highly dispersed around the galaxy. In our 
procedure for calculating the (minimum) asymmetry of each 
target we therefore shift the initial center of the asymmetry 
box by ±30 pixels in each direction, and compute individual 
asymmetries for rotation around each of the points inside 
the central 60 x 60 box. This way we obtain asymmetry 
maps of size 60 x 60 pixels for each target and filter. If 
the apparent size of the galaxy was too large to safely 
encompass the asymmetry minimum, we instead obtained 
a 120 x 120 asymmetry map. For all galaxies and all filters 
there is always a well-defined global minimum in the 
asymmetry maps, though multiple local minima exist for 
targets with highly irregular burst morphologies. In Table [8] 
we show the Petrosian minimum asymmetry in each filter 
and the corresponding Petrosian radius used to calculate 
that minimum. 

As Table [8] indicates, the Petrosian r [77(0.2)] radius is 
usually very small. It includes only the brightest isophotes, 
while faint extended irregular features like tidal tails and 
arms are thus excluded from the asymmetry measurement. 
It seems therefore prudent to have a second measure of 
asymmetry, with a different (larger) limiting radius, in 
order to capture any effect on the asymmetry that faint 
irregular structures might have. We chose the Holmberg 
radius r(fj,B = 26.5) for the B band, and in the rest of 
the optical filters the radius at which the corresponding 
surface brightness profiles reach /j ~ 26.5 mag arcsec -2 . In 



the NIR regime the H and K profiles are much shallower, 
so a limiting radius of r(fj,Nin = 23) was chosen instead. 
The area enclosed in these radii is much larger than inside 
the r[?7(0.2)] radius, however, the sample galaxies are at 
significantly differing redshifts. The difference in distance 
between the closest and the farthest galaxy is a factor of 
7. This will lower the asymmetry measured from galaxies 
further away and thus skew the asymmetry distribution 
of the sample. To even the play field and allow direct 
comparison between the galaxies, all images were first 
smoothed by a boxcar average of 1 x 1 kpc. The asymmetry 
minimum was then obtained in the same way as described 
above for the case of Petrosian asymmetries. These values 
are shown in Table [9] 



The c oncentration index was calculated following IConselicel 
(2003) as 



C = 5 x log 



rso% 

^20% 



(2) 



where r$o is the radius which encompasses 80% of the light 
inside the Petrosian 1.5 x r[T7(0.2)] radius, and 7-20 is 20% 
of rso- Measurements for each galaxy and filter are given in 
Table ITUl 



11 CHARACTERISTICS OF INDIVIDUAL 
GALAXIES 

We will adhere t o th e morphological classi fication 
of lLoose fc Thuanl l|l986l. LT86) and the ISalzer et all 
(|l989bh classification based on spectral features, and 
either reiterate the BCG c lass where available f rom 
LT86, ICairos et all |2001al ) or iGil de Paz et all (|2003l ) or, 
for the galaxies in our sample without published classi- 
fication, we will obtain a morphological classification by 
analyzing the contour and RGB plots for our data. The 
LT86 classes we apply are iE (irregular inner and elliptical 
outer isophotes), nE (central nucleus in an elliptical host), 
il, C (off-center nucleus in a cometary host), a nd iI,M (off- 
center nucleus in an apparent merger). The ISalzer et alj 
(1989b) classes, available only for UM galaxies, are DANS 
(dwarf amorphous nuclear starburst galaxies), Starburst 
nucleus galaxies, HUH (H II hotspot galaxies), DHIIH 
(dwarf H II hotspot galaxies), and SS (Sargent-Searle 
objects). 

Unless otherwise specified, the isophotal and elliptical 
surface brightness profiles agree very well. For some targets 
the isophotal profile is slightly fainter at large radii than 
the elliptical one. This is a natural consequence of the 
galaxy having a significant nebular component, but which 
is not dominating the luminosity at every radius. For 
the rest of the galaxies there is no discernible difference 
between the two profiles. This happens in the two extreme 
cases where either the dominance of the nebular emission 
is complete also at larger radii, or there is very little 
nebular emission in the galaxy. Depending on individual 
morphologies the centers of integration may be identical 
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Table 7. Estimated luminosity in excess of the exponential disk defined by h r and fio. The upper and lower 
numbers for each galaxy are for disk properties derived from the radial ranges corresponding to 24 < fj,g < 26 and 
26 < i5 28 mag arcsec - 2 , respectiv ely. Fields are left bl ank where the estimation method failed. All values have 
been corrected for Galactic extinction l lSchleeel et al,lll998l) . The % column gives the relative burst contribution to 
the total galaxy luminosity in the B band. 







% 


s* 


(U - B)„ 


(B - V)* 


(V - ft)* 


(V - J)* 


(V - If)* 


(ff - K)* 


ES0185- 


-13 


74 


15.3 ± 0.3 


-0.12 ± 0.31 


0.24 ± 0.27 


0.04 ± 0.05 


0.32 ± 0.05 


1.77 ± 0.08 


0.29 ± 0.09 






95 


15.0 ±0.1 


-0.13 ± 0.19 


0.26 ± 0.13 


0.16 ± 0.05 


0.48 ± 0.05 


1.78 ± 0.06 


0.08 ± 0.17 


ES0249- 


-31 


45 


13.1 ±0.1 




0.43 ± 0.11 


0.37 ± 0.03 


0.58 ± 0.04 


2.43 ± 0.05 


0.36 ± 0.05 






86 


12.4 ± 0.1 




0.48 ± 0.08 


0.33 ± 0.03 


0.62 ± 0.04 


2.14 ± 0.07 


0.21 ± 0.07 


ES0338- 


-04 


73 


13.9 ± 0.1 




0.37 ± 0.14 


-0.04 ± 0.01 


-0.74 ± 0.05 


0.74 ± 0.03 


0.31 ± 0.04 






94 


13.6 ± 0.2 




0.35 ± 0.22 


0.05 ± 0.04 


-0.34 ± 0.06 


0.93 ± 0.05 


0.07 ± 0.13 


ESO400- 


-43 


76 1 


14.6 ± 0.2t 








0.31 ± 0.22 


1.83 ± 0.23 


0.35 ± 0.08 






94t 


14.6 ± 0.3t 








0.35 ± 0.33 


1.88 ± 0.33 


0.27 ± 0.28 



ESO421-02 



ESO462-20 


80 
95 


14.7 ± 0.2 
14.5 ± 0.3 


-0.07 ± 0.23 
-0.09 ± 0.36 


0.22 ± 0.17 
0.23 ± 0.33 






2 
2 


09 ± 
00 ± 


20 
22 


0.56 ± 
0.44 ± 


21 
27 


HE2-10 




53 
89 


12.6 ± 0.1 
12.1 ± 0.1 


0.04 ± 0.15 
0.13 ± 0.18 


0.22 ± 0.07 
0.44 ± 0.13 






2 
2 


98 ± 
78 ± 


06 
06 


0.51 ± 
0.45 ± 


09 
09 


HL293B 




14 
16 


19.5 ±0.1 
19.4 ±0.2 


-0.76 ± 0.16 
-0.60 ± 0.27 


0.42 ± 0.07 
0.05 ± 0.21 


-1.32 ± 0.07 














IIZW40 




68 
87 


12.1 ± 0.2 
11.8 ±0.9 


-0.56 ± 0.97 


0.09 ± 0.20 
0.23 ± 0.95 




-0.10 ± 0.09 
-0.16 ± 0.10 



1 


78 ± 
08 ± 


07 
30 


-0.16 ± 
0.43 ± 


19 
30 


MK600 




49 
76 


15.9 ±0.1 
15.4 ±0.1 


-0.46 ± 0.20 
-0.40 ± 0.18 


0.35 ± 0.14 
0.41 ± 0.10 


0.01 ± 0.07 
0.25 ± 0.07 


0.18 ± 0.10 
0.44 ± 0.10 


1 
1 


45 ± 
77 ± 


06 
11 


0.52 ± 
0.57 ± 


09 
16 


MK900 




49 
73 


15.0 ± 0.1 
14.6 ±0.1 


-0.37 ± 0.16 
-0.23 ± 0.17 


0.42 ± 0.07 
0.54 ± 0.10 


0.35 ± 0.07 
0.43 ± 0.07 


0.60 ± 0.09 
0.71 ± 0.09 


1 
1 


95 ± 
88 ± 


11 
12 


0.06 ± 
-0.08 ± 


10 
12 


MK930 




52 
94 


15.7 ± 0.1 
15.1 ± 0.2 


-0.80 ± 0.19 
-0.78 ± 0.22 


0.49 ± 0.13 
0.42 ± 0.17 


-0.25 ± 0.07 
0.05 ± 0.08 


-0.33 ± 0.10 
0.14 ± 0.11 



1 


94 ± 
81 ± 


10 
31 


1.17 ± 
0.75 ± 


13 
35 


MK996 




-12 
68 


16.0 ±0.1 
15.5 ±0.1 


-0.54 ± 0.15 
-0.22 ± 0.21 


0.45 ± 0.05 
0.54 ± 0.15 


0.33 ± 0.06 
0.40 ± 0.06 


0.49 ± 0.09 
0.70 ± 0.09 


2.54 ± 0.05 
2.44 ± 0.09 


0.85 ± 


09 


SBS0335- 


052E 


75 


17.1 ±0.1 








-0.14 ± 0.11* 


0.75 ± 0.24* 


0.50 ± 


24 






86 


16.9 ±0.3 








-0.19 ± 0.32* 


OA 


9 ± 0.37* 


0.44 ± 


26 


SBS0335- 


052W 


26 


20.6 ± 0.1 








-0.11 ± 0.10* 


0.57 ± 0.22* 










39 


20.2 ± 0.3 








-2.33 ± 0.29* 


0.5 


2 ± 0.56* 






TOL0341 


-407 


83 


16.6 ±0.1 


-0.50 ± 0.19 


0.65 ± 0.12 




0.01 ± 0.50 


1 


40 ± 


08 


0.25 ± 


10 


TOL1457 


-262 I 


73 
94 


14.5 ±0.2 
14.3 ± 0.2 




0.29 ± 0.18 
0.27 ± 0.22 






1 
1 


26 ± 
22 ± 


11 
10 






TOL1457 


-262 II 


47 
57 


15.6 ±0.1 
15.9 ± 0.5 




-0.06 ± 0.10 
-0.65 ± 0.45 






1 


91 ± 


10 






UM133 




56 


16.1 ±0.1 


-0.61 ± 0.18 


0.26 ± 0.10 


0.18 ± 0.06 


0.09 ± 0.09 












UM160 


UM238 




21 

60 


18.6 ±0.1 
17.5 ± 0.2 


-0.40 ± 0.18 
-0.32 ± 0.24 


0.38 ± 0.11 
0.34 ± 0.19 


0.31 ± 0.07 
0.44 ± 0.07 


0.83 ± 0.10 
0.88 ± 0.10 


1 


97 ± 


08 


0.77 ± 


15 


UM417 




23 


19.8 ±0.1 


-0.85 ± 0.15 


0.80 ± 0.16 


-0.70 ± 0.09 


-1.55 ± 0.09 


-0 


29 ± 


17 


0.71 ± 


25 


UM448 




78 
88 


14.4 ± 0.2 
14.3 ±0.1 




0.47 ± 0.23 
0.46 ± 0.10 


0.20 ± 0.08 
0.24 ± 0.06 


0.53 ± 0.07 
0.59 ± 0.06 


2.43 ± 
2.50 ± 


13 
20 






UM619 




60 


16.3 ± 0.3 


0.08 ± 0.28 


0.20 ± 0.25 


-0.36 ± 0.07 


0.21 ± 0.07 





77 ± 


10 


-0.22 ± 


17 



f - V instead of B filter 
% - B instead of V filter 



or vary in physical location, which we will mention in 
the cases where it has any bearing on the analysis. We 
interp ret the color radial profiles using the iMarigo et al.l 
(|2008f > stellar evolutionary tracks for a Salpeter IMF stellar 
population with metallicities Z = 0.001, 0.004, 0.02, 0.2. 
These tracks are predominantly used when comparing with 
the old host population. Here we will briefly mention the 
main guidelines such evolutionary tracks can offer. Some 
combinations of colors, e.g. B — V vs. V — I, suffer from 
a severe age-metallicity degeneracy, getting simultaneously 
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redder with increased age and with increased metallicity 
in a fashion which often makes it impossible to distinguish 
between the two. For old populations (> 1 Gyr) combina- 
tions involving V — K are useful in breaking the metallicity 
degeneracy because V — K is very weakly dependent on 
age but has a strong Z dependence, getting redder with 
increased Z. For all color combinations there are limiting 
cases which can sometimes be useful to indicate Z and age 
with a fair certainty, e.g. V — I > 1.1 always indicates a 
high Z (close to solar or above) and an age older than 5 
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Table 6. Host structural parameters. Absolute B magnitude Since BCGs are anticipated to have a significant nebular 

of the host (M B ost ) obtained by subtracting the burst luminosity. component, we also look at stellar evolutionary tracks 

The scale length h r in arcseconds and kpc, and the central surface with added nebular emission. These mod els are based on 

brightness [mq are based on a weighted least squares fit to the the Yggdrasil spectral synthesis code (jZackrisson et all 

deepest image (B band) for two radial ranges derived from 24 < |201lf) w ith Starburst99 Pad o ya-AGB stellar population 

fi B < 26 a nd 26 < n B < 28 mag arcsec" 2 . Correction for Galactic spectra (jLeitherer et all 1 19991 : IVazauez fc Leithererl l2005h 



extinction J Schlegel et al. |(l99i) has been applied. at metallicities Z = . 004, Z = 0.004, Z = 0.008 and 

z = 0.020. The nebular contribution to the overall spectral 

Galaxy >'b Ajj|° st \ >v pc e nergy distrib ution is computed u sing the photoionization 

Esoi85-i3 24-26 -i7.7 4.18 ±o.26 1.55 ± o.io 20.87 ±o.27 code Cloudy l|Ferland et al.l 11998). assuming a spherical 

-geometry for the photoionized gas and a nebular metallicity 
identical to that of the stars. With the help of these two 

ESO33S-04 ~ ~ 7.74 ±0.22 i.4o ±0.04 20.68 ±o.i4 models we carr y out a separate analysis for each galaxy 
26-28 -i6.2 18.93 ±1.36 3.43 ±0.25 23 .84 ± o. 22 anc j tne i r respective total, central (down to fi B ~ 24 



26-28 -15.9 15.97 ± 0.94 5.91 ± 0.35 24.78 ±0.12 

ES0249-31 24-26 -17.0 31.48 ± 0.87 1.65 ± 0.05 21.40 ±0.11 

26-28 -15.7 74.85 ± 1.92 3.92 ± 0.10 24.23 ±0.08 



Eso4oo-43t 24-26 -i8.4 4.3i ±0.10 i.65 ±0.04 20.64 ±0.22 mag arcsec ), and burst colors, as well as the colors 

26 - 28 - 16 ° 3-39 ± 0.75 3.61 ± 0.29 23.60 ±0.32 p f ^ ^ = ^ _ ^ ^ = ^ _ 2 g mag ^^-2 

host regions. The total number of separate color-color 
-diagrams we have analyzed is ~ 450, and we do not 



ESO421-02 24-26 11.97 ± 0.09 0.72 ± 0.01 20.57 ± 0.05 

26-28 11.11 ± 0.28 0.67 ± 0.02 20.20 ±0.18 



ESO 162-20 24-26 -17.9 6.33 ± 0.30 2.44 ± 0.12 21.27 ±0.17 



26-28 -16.7 16.49 ±2.33 6.36 ±0.90 24.29 ±0.33 present those due to space considerations, but the plots are 

"available on demand. From the individual analysis of the 
_color-color diagrams we will attempt to give indications of 



24-26 -18.0 17.94 ± 0.22 1.37 ± 0.02 20.49 ± 0.07 

26-28 -16.6 41.13 ± 1.59 3.15 ± 0.12 23.53 ±0.13 



24-26 -12.7 2.25 ± 0.01 0.18 ±0.00 

26-28 -12.9 2.28 ± 0.08 0.18 ± 0.01 20.73 ± 0.20 



20.66 ± o.o5 the age and metallicity of both the underlying host and 
-the young populations for each galaxy. An example color- 



24-26 -17.1 29.26 ± 1.83 1.67 ± 0.10 21.60 ±0.20 . .. , . rrrn r ,i i i , 

26-28 color diagram can be seen m b lgureyjjior the whole sample. 



24-26 -13.9 11.56 ± 0.41 0.61 ± 0.02 21.58 ±0.12 

26-28 -13.4 25.94 ± 0.70 1.37 ± 0.04 23.79 ±0.09 



We have briefly looked at Ha equivalent width radial 



'el's ~Atl ItlZia.ll i.Toia.ol 23:03 lo.ol P rofiles for ei S ht of the galaxies in the sample, however, 
■ detailed analysis of Ha data is beyond the scope of this 

24-26 -18.3 3.84 ± 0.07 1.33 ± 0.02 19.59 ±0.11 

26-28 -16.0 11.30 ±0.61 3.9i ±0.21 23.85 ±0.16 paper. Here we only note that if nebular emission is 
24_26 _i5.4 5.91 ± 0.03 o.52 ± 0.00 20.69 ±0.03 detected at large radii, then there is a tendency for the 

26-28 -14.9 7.35 ± 0.20 0.65 ± 0.02 21.78 ±0.14 



_equivalent width to be significantly large (> 400-~ 2000A) 



SBS0335-052E 24-26 -i4.9 2.59 ±0.04 0.68 ±0.01 21.66 ±0.07 at a u rac iii meaning the colors in every radial bin we 

26-28 -14.7 3.93 ± 0.30 1.03 ± 0.08 23.07 ± 0.30 ' J 



-study will have nebular emission contribution. For galaxies 



SBS0335-052W 24-26 -13.2 2.62 ± 0.01 0.69 ± 0.00 22.60 ±0.05 . , , , r ... , , 

26-28 -13.8 2.81 ± 0.20 0.74 ±0.05 22.95 ±0.27 with central star forming regions the ettect or the nebular 

2.54 ±0.13 0.75 ±0.04 18.71 ±0.32 component is negligible at large radii (^ B = 26-28 

26-28 -i5.4 7.79 ±0.26 2.29 ± 0.08 23.56 ±0.11 m aR arcsec -2) if ne bular emission is indicated by the 

24 - 2 6 5.90 ±0.25 3.63 ± o.o7 20.84 ±0.18 stellar evolutionary models even at these radii, the Ha 

26-28 -18.3 12.55 ± 0.87 3.39 ± 0.22 23.62 ± 0.22 J ' 

-data show that it is in-situ for the galaxies we have looked at. 



TOL0341-407 24-26 



TOL1457-262 / 



TOL1457-262// 24-26 -18.5 10.29 ± 0.20 2.08 ± 0.09 20.46 ±0.10 

26-28 -16.9 9.61 ± 0.63 4.43 ± 0.31 20.08 ± 0.45 



24-26 9.84 ± 0.25 0.87 ± 0.02 20.89 ±0.11 

26-28 -14.8 15.53 ± 0.33 1.38 ± 0.03 22.80 ±0.10 



In our data the color profiles for most targets show a radial 
_reddening which is likely to be mainly due to the increase of 



24 ~ 26 10.36 ±o.62 i.4i ±0.08 20.27 ±0.28 stellar age at large radii, unless otherwise specifically men- 

26-28 7.80 ± 0.28 1.06 ± 0.04 18.70 ± 0.30 ° ° ' 



24-26 -16.0 5.20 ± 0.14 1.37 ± 0.04 21.51 ± 0.10 

26-28 -15.' 



-tioned. In some cases we observe a flattening of the color 



8.30 ± 40 2 i8 ± a.ii 2311 ± 0.18 profiles beyond some radius. This is expected in the case of 

a sustained low star-formation activity randomly percolat- 



24-26 -12.7 7.28 ± 0.54 1.19 ± 0.09 22.96 ±0.14 

26-28 3.10 ±0.12 0.51 ± 0.02 19.76 ±0.26 



26-28 -17.9 23.18 ± 0.89 9.27 ± 0.36 24.33 ± 0.09 

24-26 4.22 ± 0.11 1.40 ±0.04 19.18 ±0.16 

26-28 -17.4 6.45 ± 0.29 2.14 ± 0.10 21.48 ±0.25 



f - V instead of B filtc 



ing across the disk w hen there are no favored locations in 

24-26 -i8.i 10.73 ±0.80 4.29 ±0.32 22.53 ±0.22 terms of gas density (van Zee 200l]). In other words, such 

profiles indicate that we are sampling a population which 
is homogeneous in the sense that its different constituents 
%ave similar metallicities and ages. Note that this assumed 
homogeneity applies to the host, and not to the composite 
galaxy. The flattening of the color profiles is more easy to see 
in optical-NIR colors than in the purely optical, for which 
Gyr. Similar limits exist for the rest. H — K is the only the host is usually sampled at much larger radii. In some 
color that beyond 1 Gyr gets bluer with increased age for cases the isophotal color profiles are better tracers of radial 
any Z but the tendency is weaker for solar and super-solar age or metallicity evolution than the elliptical integration 
Z. ones because the irregular morphology of the targets will di- 

lute any accumulative effect that the e.g. increasing age has 
In all cases nebular emission and the p r esenc e of dust on the profile shape. The effects of extinction occasionally 
modify the observed colors, e.g. iHunt et al.l (|2003l ) find that show up in the region of the burst, typically reddening the 
H — K can redden to ~ 0.5, while V — I can get ~ 0.8 innermost parts of the profile at and around the brightest 
mag bluer if the region is dominated by nebular emission. isophotes. In the event that this is not observed, compact 
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starburst knots would typically still give rise to strong fluc- 
tuations in the colors at small radii and those colors will 
be strongly affected by nebular emission. We will be more 
prolific in our analysis for galaxies which have been sparsely 
studied in the literature, or for which there is little or no 
analysis of the radial surface brightness and color distribu- 
tions. 



ES0185-13 

This strongly starbursting galaxy (|Ostlin et al.l l200ll . 
SFR = 5.4 Mq/jx) is likely the product of collision as 
evidenced by the detection of numerous bright compact 
star clusters (|6stlinl l2000bh and by the galaxy morphol- 
ogy. Beyond /is ~ 26.0 mag arcsec -2 the elliptical surface 
brightness profile reveals a prominent low surface bright- 
ness (LSB) component, which was previously undetected. 
In Appendix IB 1 1 we show how much detailed structure of 
this component we actually detect. Clearly visible in all op- 
tical filters is the flux overdensity to the North-East of the 
prominent tail. This is not seen in the NIR data, which is 
not deep enough to reach such large radii. In the NIR we 
only detect a single exponential disk (with no central ex- 
cess) down to fiH,K = 23 arcsec -2 at r ~ 15 arcsec, which 
is where the B band profile break occurs. We classify this 
galaxy as iI,M as it is an apparent merger. The color-color 
diagrams indicate the presence of an old population (> 2 
Gyr), as well as strong nebular emission contamination for 
the star-forming regions, which drops considerably beyond 
(iB > 24 and becomes negligible for fj,B > 26 mag arcsec -2 . 
While we were unable to constrain the age of the starburst 
in this case, it is quite young, as i ndicated by the 3.5 Myr 
old peak of star cluster formation (|Adamo et al.ll201ll 1. 



as iE BCG, but its chaotic velocity field is signi ficant indica- 
tion that this galaxy is in fact a recent merger (|6stlin et al.l 
I2001T) . It has a spectrosc opically confi rmed companion at 70 
kpc projected distance l|6stlin et al.lll999h . Strong nebular 
emission contamination is visible for all optical colors, but 
the effect is negligible for the old population region (p,s > 26 
mag arcsec -2 ). Optical-NIR and pure NIR colors detect the 
presence of a young < 100 Myr stellar population. Note 
that there is a foreground star South-West of the central 
star-forming region, which we have masked out before any 
integrated colors or profiles were obtained. 



ESO400-G043 

There are two prominent star-forming regions in this galaxy, 
both of which are off-center. We classify it as iE BCG due to 
its regular elliptical outer isopho tes. It has a dl companion 
at a projected distance of 70 kp c jBergvall fc Jorsaterl 19881) 
with a co mparable H I mass. iBergvall fe Olofssonl ( 1984 . 
1986) and IBergvall fc Ostlinl (|2002T ) demonstrate the pres- 
ence of an old stellar population. We cannot completely 
isolate this population in any of the investigated regions, 
though the measured V — I > 0.5 at large radii would re- 
quire stars > 1 Gyr of age. Additionally, the strong negative 
H — K gradient we observe indicates a radial decrease in 
metallicity and/or an increase in age. The color-color di- 
agrams indicate that the integrated colors for all regions 
contain young and old stars, mixed with gas. Specifically, 
contribution from gas is clearly visible in V — K which is 
too blue to be com ing solely from stars. Th e dust content in 
this galaxy is low (|Bergvall fc Ostlinl [2002] '). 



ES0249-31 

This peculiar spiral with ex tended arms shows clear signs 
of an ongoing merger event |Bergvall et al.ll2003T l. We thus 
classify it as iI,M BCG. The profiles should be affected 
by a significant amount of extinction (jBergvall et al.l 120031 . 
H a /Hp ~ 4.1). The star-forming regions are spread out at 
various radii resulting in many small scale fluctuations in 
the optical color profiles but no clear large-scale gradients. 
In hybrid and purely NIR colors, both the radial color pro- 
files and the integrated colors in color-color diagrams show 
a negative radial metallicity gradient. The presence of an old 
stellar population > 2 Gyr with low or intermediate metal- 
licity (Z < 0.008) is indicated, as well as a younger popula- 
tion with metallicity closer to solar (0.008 < Z < 0.02). 



ESO338-I04 



This Wolf-Rayet galaxy jSchaerer et all 1 19991 ) has been 
thoroughly stu died in the optical broa d band and nar- 
rowb a nd (e.g. IBergvall fc Olofssonl 1 1984 1 19861 ; iBergvalll 
1 19851 ; IBergvall fc Ostlinl |2002| V It has undergone previ- 
ous and rather recent star-formation events, as evidenced 
by i ts rich population of intermediate-age globular clus- 
ters |6stlin et al.lll99"8h . which are unresolved in o ur images. 
Based on its morphology, iGil de Paz et all (|2003l ) classify it 



ESO421-02 

To gether with ESQ462-IG 02 0, this galaxy was first studied 
by Bergvall fc Olofssar] l[l986'). It shows a spatially extended 
starburst with multiple nuclei, embedded in a common el- 
liptical host. We therefore classify it as iE BCG. Neither 
of our population synthesis models fits the observed colors 
for the burst or the two outer regions. This is probably due 
to contribution from nebular emission which is not in-situ, 
which the models do not account for. The color profiles are 
predominantly flat. U — B shows gradients, but this color 
is very sensitive to an aging young population and would 
reflect such aging by changing on shorter time scales than 
the redder colors. Flat color profiles, minor local color vari- 
ations notwithstanding, are a signature of a percolating SF 
at a low but persistent rate (van Zee 2001). However, the 
RGB image suggests that this is a merger since two very 
similar-looking dwarf galaxies are seen in physical contact 
with each other, with the starburst located between them 
along the region of contact. This is further supported by 
unpublished kinematic data, which give a clear indication 
that this is indeed an on-going merger between two dwarfs, 
one to the East, the other to the West. The Eastern part is 
falling into the Western part because the too regular (blue) 
velocity curve of the latter, implies that it could not have 
suffered a fly-through of the Eastern dwarf (which has a red 
velocity curve). 
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ES0462-IG 020 

This galaxy has been classified as a distorted contact pair 
in the surface photometry catalog of ESO-Uppsala galaxies. 
We see indications of the ongoing merger in the optical and 
optical-NIR isophotal color profiles, which show sharp color 
changes visible as discontinuity jumps at r ~ 1,2, and 4 
kpc. These jumps in color imply the presence of more than 
two distinct stellar populations with different SFH which 
have mixed spatially on large but not on small scales. We 
therefore classify this galaxy as iI,M BCG. The color-color 
diagrams indicate that perhaps more than two populations 
are present and influencing the integrated colors. Nebular 
emission is implied, but the current SF is poorly modeled 
with an instantaneous burst. The only definite conclusion 
we can draw from the color-color diagrams for this galaxy is 
that a population older than 1 Gyr is present, but we cannot 
narrow down its metallicity. There is a bright nearby star to 
the North-East of the galaxy which could be contaminating 
the galaxy light beyond /ib ~ 26 mag arcsec -2 . 



HE2-10 (ESO495-G021) 

He2- 10 is the first Wolf-Rayet galaxy ever discov- 
ered (|Allen et al . 1976). Contrary to popular belief this is 
not a dwarf as it is too luminous {Mb = —19.1, Mk = 
—22.3 mag) and has an extended regular envelope spanning 
~ 13.8 kpc. The observed V — K color profile is centrally 
very red, r eflecting both a hi g her Z , measured to be close 
to solar by iKobulnickv et all (Il999l ) , and a larger internal 
extinction. There is dust clearly visible in the RGB image 
(dark lanes in the North and South directions). The negative 
V—K gradient is thus due to a strong metallicity and extinc- 
tion decrease with radius. In fact we can model the regular 
envelop surrounding the starburst with a metal-poor stellar 
population of Z ~ 0.004 - a sharp decrease from Zq. The 
H — K profile on the other hand is essentially flat beyond 
~ 15 arcsec. It is also too red at any radius for any age 
and metallicity, as is ofte n the case for Wolf-Rayet galax- 
ies (|Buckalew et alj|2005l ). Its mean value of 0.5 requires a 
large contribution from nebular continuum. If we are sam- 
pling a population older than ~ 4.5 Gyr beyond the central 
1.5 kpc then the expected radial decrease in Z, gas density, 
and dust content, together with the simultaneous increase 
in age will all enhance each other to produce a strong neg- 
ative gradient. The flatness we are seeing must therefore be 
due to an old homogeneous population, or to a population 
of an intermediate or younger age, which is well-mixed. The 
regular contours of the outer isophotes support the former 
case. We classify this galaxy as iE BCG 



HL 293B (SDSS J223036. 79-000636. 9) 

This faint dwarf is a Haro-Luyten object, quasi- 
stellar in appe arance. It was first classified as a BCG 



spatiall y coextensive with the burst nucleus. ICairos et alj 
l)2001ar ) similarly do not find an extended underlying compo- 
nent in this galaxy. There is nevertheless an old stellar pop- 
ulation present, although it does not appear to be spatially 
dominant. While all other colors show signs of strong con- 
tamination by nebular emission, V — K > 2.0 at V — I > 0.5 
would require the presence of old stars. 



IIZw40 

This iI,M BCG galaxy is very blue and strongly starburst- 
ing ()Vanzi et al] I1996T ). Its op tical and NIR spectra have 
been extensively studied (e.g. iBaldwin et al] 1 19821 ; iThuanl 
1983h. as well as its m o rphol ogy and companions (e.g. 



by iThuan fc Martini dl98lf>. which was later refined to nE 
BCG bv lCairos et al] 1 2001a! '). This galaxy is currently un- 
dergoing a starburst only a few Myr old, as evidenced by the 
detec tion of luminous blue variable stars (jlzotov fc Thuanl 
2009). It seems to have a lot less NIR than optical flux, even 
though the exposure times are reasonably large (42 minutes 
in K). The NIR flux is coming from a very compact region 



Campos- Aguilar fc Moles! Il99lf ). H I gas content (e.g. 
■ Sage et al] 1992 ) and met a llicity (e.g. Kunth fc Sargent] 
1 19831 : iTerlevich et all Il99ll : iMarconi et al.l Il994l ). IIZw40 
is the result of a merger of two dwar f galaxies, judg- 
ing by its optical morpho logy (e.g. ISargent fc Searlel 
1 19701 : iTelles fc Terlevichl 1 19971 ) and its extreme H I tidal 
tail extending South-E ast approximately 8.5' in diame- 
ter (|van Zee et al.| [l998) which does not seem to have an 
optical counterpart. The central colors are consistent with 
a very young stellar population of a few Myr, while at 
large radii all colors are strongly contaminated by nebular 
emissi on. The dust mass found in this galaxy is rather 
large dLisenfeld fc Ferraral [i~99Sl . M du3t » 10 6 M W ) which 
is comparable to the stellar mass (e.g. IVanzi et al] Il996l . 
M* « 2.5 x 1O 6 M ), so the reddening of the isophotal 
U — B in the starburst regions is likely the effect of internal 
extinction. 

There is no radial range in our data at which we can sam- 
ple a region dominated by an old population. Since IIZw40 
has such low Galactic latitude it suffers great extinction (3.5 
mag in B) and the otherwise quite long exposure times of e.g. 
40 minutes in the B band only allow us to reach \ib ~ 25.0 
mag arcsec -2 . Judging from the rest of our sample, the sur- 
face brightness profile in the region around [is ~ 25 mag 
arcsec -2 is still very much contaminated by the burst, espe- 
cially if the galaxy has disturbed morphology. It is therefore 
entirely possible that we are only seeing the current young 
starburst and, at large radii, the well-mixed remnants of 
some intermediate-age population from a previous burst. 
There could be an old host dominant beyond /j,b ~ 26 mag 
arcsec -2 , as is the case for e.g. ES0185-13 and UM448, and 
this would be entirely invisible to us. 



MK600 

At first glance this iE BCG has a double exponential disk 
profile, with a faint population which dominates the sur- 
face brightness profile beyond /is ~ 26 mag arcsec -2 . Pre- 
vious studie^_oi]_tius_gal^y_do not mention such a struc- 
ture l|Hunter fc Elmegreenll200r3 '). Our B band profile is the 
deepest to date, however, MK600 is situated between three 
very bright stars and the second exponential disk compo- 
nent could be due to contamination from these stars. This 
does indeed seem to be the case judging by the deep contour 
image in Appendix IB3I We have applied heavy masking of 
these stars. The profile data beyond fiB ~ 26 mag arcsec -2 
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are based on a narrow unmasked cone to the South- West, 
with minimum cone radius ~ 1' at fis ~ 26 mag arcsec -2 . 
As this area is seemingly unaffected by the stars, this struc- 
ture could still be real, but we remain suspicious of its exis- 
tence and will not dwell on it further. In the H and K bands 
we see a single exponential disk profile down to (J,h,k = 23 
mag arcsec -2 at radius r ~ 25 arcsec, with no indication of 
a central excess. lAugarde et alj (|l990T ) find strong evidence 
against a single starburst episode in MK600, and favor in- 
stead a mode of sustained star formation over ~ 100 Myr. 
We find some supporting evidence for this since the color- 
color diagrams strongly indicate that the current star for- 
mation is poorly modeled with an instantaneous burst. The 
galaxy is a faint IRAS source dBicav et ail 1 19951 ). with low 
extinction (|Hunter fc Elmegreenl |2004 E(B - V) = 0.05), 
and indeed there does not appear to be a lot of dust ob- 
scuring the young blue regions, as evidenced by the location 
of the burst colors in a U — B or B — V vs. H — K dia- 
grams, which are both consistent with negligible amounts of 
dust. Nebular emission contribution is visible when compar- 
ing B — V with U - B, V - R, V - I, or V - K at all radii, 
though the effect is considerably smaller beyond 20 arcsec. 



MK900 

The star forming regions in this nE BCG seem to follow 
the semi-minor axis of the underlying galaxy and are found 
symmetrically on both sides of the geometric center. All 
color-color diagrams ubiquitously indicate the presence of 
a homogeneous old (> 4 Gyr) metal-poor (Z < 0.004) host. 
In this case the instantaneus burst assumption fits the ob- 
served colors very well. U—B shows no extinction even in the 
SF regions, w hich is to be expe cted as this is another faint 
IRAS source l|Bicav et ail 1 19951 ). Contamination from neb- 
ular emission is confined to the star formation region, with 
its effect on the colors being negligible beyond the central 
10 arcsec. Thi s is consistent w ith the galaxy being moder- 
ately gas-rich (|Wei et al.ll2O10l . M H i +He/M i , = 0.34). The 
nuclear star formati on has an H a luminosity comparable to 
that of 30 Doradus (|van Zeell2o6ch . 



MK930 



This is a Wolf-Rave t galaxy JSchaerer et al] 1 19991 ) whose 
multiple nuclei (e.g. iGimeno et al.ll2004h ' indicate the ad- 
vanced stages of a merger. We therefore classify it as an 
iI,M BCG. The presence of an old (> 4 Gyr) stellar popu- 
lation can be deduced from all color profiles. Nebular emis- 
sion is implied by all color-color diagrams, but at the same 
time neither of the models is well-suited for the burst, to- 
tal, or central mixed colors, indicating either a prolonged 
burst, or severe dust effects. The latter seems more likely 
since non-ne gligible amounts of warm dust are present in 
this galaxy (|Wu et al.l 120081 ). We were unable to impose 
usefu l constraints on the a ge of the young population, how- 
ever. |Xd amo et all (|201 lh place it at < 10 Myr based on 
the ages of super star clusters. Note that the flattening of 
the elliptical B band profile beyond fj,B ~ 27 mag arcsec -2 
is due to the sky, and is not an indication of an additional 
underlying component. 



MK996 

This nE BCG is also a Wolf-Rayet galaxy with a com- 
pact off-ce nter star formation region. It is old, as evi- 
denced by iThuan et ail l|l996l ) who find an old globular 
system (~ 10 Gyr) distributed around the galaxy. All of 
our data are consistent with an old (> 3 Gyr) metal-poor 
(Z < 0.004) stellar population with moderate contamina- 
tion of nebular emission predominantly visible in U — B, 
B — V, and H — K diagrams, which is negligible outside the 
compact star formation region. The pr esence of a du s t lane 
has previously been confir med by both lThuan et all (| 19961 ) 
and boublier et al l (| 19991 ). and is also visible in our RGB 
image. 



SBS 0335-052 E 

The K band exposure time for this galaxy is the longest 
one in the entire sample (~ 5 hours with NOTCAM), yet 
we detect very little K band flux. The surface brightness 
quickly drops to fix ~ 22.0 mag arcsec -2 within the central 
2 kpc (r « 4 arcsec), beyond which radius the errors become 
too large (a > 1 mag). We performed a sanity check using 
the ~ 2 hour SOFI K band image, but the result was the 
same and we could only reach fijf ~ 22.0 mag arcsec -2 . The 
K band flux is concentrated in the compact regular ellipti- 
cal region visible in the contour plot, with no measurable 
flux beyond this region. Considering this it is difficult to 
speak of a host galaxy in the traditional sense of the term 
since the B band, which contains contributions from both 
stars and gas, seems to be much more extended than the 
K band, which is considered a good proxy for stellar mass. 
The radius r ~ 2.5 arcsec marks the "end" of the compact 
K band flux region. The abrupt onset of this change testifies 
that the two regions contain populations with very different 
physical properties. A very strong contribution by nebular 
emission is necessary to reconcile the observed H — K colors 
with B — I or B — K, the latter two colors indicating a very 
young population of a few 10 Myr. This is consi stent with the 
~ 10 Myr age of super star clusters found bv lAdamo et al.l 
(|20ld) . Indeed, rese archers have argued that this is a pri- 
mordial galaxy (e.g. 



iPapaderos et a l. 1998; Vanzi et al 



Izotov et al 



1997 



IThuan et alj Il997l ; 
2000) with the estimated 



age of the exten ded component dra matically decreasing from 
~ 100 M yr fe.g.llzotoy et al.lll997l ) to 5-6 Myr in more recent 
resea rch l|Vanzi et al. 2000l ) (see, however. lOstlin fc Kunthl 
l200ll for a contrary o pinion). It also has surprisingly large 
amounts of hot dust l|Thuan et ail [19991 ; I Vanzi et al.l |200 ~ 
3-5 x 1O 5 M ). 



SBS 0335-052 W 

This galaxy was first disco vered as a peak i n the H I maps 
of the companion galaxy. iPustil'nik et al.l (1 19971) detected 



(|2005ls measured 



an optical counterpart, and llzotov et al 
its record-low metallicity, ~ 10% lower than the previous 
record-holder I Zw 18. This galaxy is in the same field of 
view as SBS 0335-052 E, so the same remarks on K band 
exposure time apply here. Similarly, we find very little K 
band flux, again located only in the small compact region 
visible in the contour plot. The K band profile reaches \ik ~ 
22.5 mag arcsec -2 but we can barely resolve the galaxy in 
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this filter - there are only 3 reliable data points in the profile. 
Due to this we shall not attempt to interpret the B — K and 
H — K color profiles. The B — I profile behaves similarly 
to SBS 0335-052 E, with a pronounced blue minimum at 
r ~ 1.5 arcsec, which, again, marks the "end" of the compact 
bright region seen in the contour plot. 

Tol 0341-407 

The blue star-forming regions in the East and West are two 
separate galaxies wit h distinctly differen t photometric and 
kinematic properties (|Qstlin et al.lll999h . hence we classify 
this merger as iI,M. Due to the very irregular morphology of 
the galaxy the isophotal and elliptical centers of integration 
are very different and the two surface brightness profiles do 
not trace each other at any radius. The elliptical color pro- 
files start at the center of mass (using the center of light 
as proxy), which happens to be close to the geometric cen- 
ter due to the axisymmetric placement of the two starburst 
regions. Judging hyB — VoxU — BvsV — I diagrams, 
the central region is completely dominated by a population 
younger than 1 Gyr, while beyond (j,b ~ 26 mag arcsec -2 
we detect another population, > 5 Gyr. Thus the color- 
color diagrams are consistent with the existence of three 
separate stellar populations in this galaxy, as already sug- 
gested bv lWestera et all (|2004l ) . Strong nebular emission is 
contaminating the colors of the two starbursting regions in 
U — B vs B — V and B — V vs V — K diagrams. 

Tol 1457-262 J 

This is a merger between two distinct galaxies, Toll457- 
262a (North) and Toll457-262b (South). We thus classify 
it as iI,M. The lack of U, I and H filters prevents us from 
imposing any kind of limit on gas content, stellar age or 
metallicity. The burst region is spatially dominating and 
shows strong nebular emission contamination in a B — V 
vs. V — K diagram. The /ib = 24-26 mag arcsec -2 region 
shows an older metal-rich (~ Zq) population, however we 
cannot confirm this with only three filters. Note that here 
the K band is suspect in the faint regions, since upon heavily 
smoothing the K band data to closer examine the sky we de- 
tected a residual background structure, which is not additive 
in nature, and clearly indicates an illumination/flatfielding 
problem. This residual structure completely dominates the 
K band profile beyond ~ 18 arcsec for Tol 1457-262 I and 
beyond ~ 35 arcsec for Tol 1457-262 II. 

Tol 1457-262 77 

At first glance it looks deceptively easy to fit this object 
with elliptical isophotes, however, the isophotal and ellipti- 
cal surface brightness profiles are very different from each 
other, though with similar centers of integration. The rea- 
son for this discrepancy becomes obvious once we compare 
the distances to the different components of this object. 
Toll457-262c (roughly, the Eastern region) and Toll457- 
262e (the Western region) have a tabulated NED velocity 
difference of ~ 600 km s -1 . In other words, these are two 
different galaxies superposed along the line of sight with 
a velocity different large enough to make it unlikely that 



they are interacting with each other since this is signifi- 
cantly high er than the li mit for Hickson compact groups 
of galaxies (|Hicksonlll997h . The galaxies are superposed in 
a fashion which makes it impossible for us to separate them 
using just broadband imaging data. As a consequence, both 
surface brightness profiles are sampling a mixture of stellar 
populations at different redshifts. 

UM133 

This is an U,C BC G, with no evidence of interac- 
tion jTelles et al.lll997l ). The elliptical color profiles are in- 
fluen ced by the off-center location of the strongest starburst 
knot. lNoeske et alj <|2000h suggest that cometary BCGs like 
UM133 and UM417 could be an evolutionary link between 
BCGs and other dwarf galaxies because their structural pa- 
rameters (no and Hr) are intermediate between those in- 
ferred for iE/nE BCGs and the dlrrs/dEs. Color-color dia- 
grams indicate moderate contamination from nebular emis- 
sion. Due to the morphology of the galaxy we cannot give 
separate metallicity and age limits for the different popula- 
tions, but we do detect at least two, one > 2 Gyr, and the 
other younger. 

UM160 

Numerous very blue star forming regions are spread out 
aroun d the red center of this Wolf-Rayet l|Schaerer et al.l 
1999) galaxy. The starbursting knots are tracing what looks 
like distorted spiral arms. This is a gas- rich (~ 9 x 1O 9 M ) 
galaxy with a low-mass ( ~ 3 x 10 8 Mp)) dwarf comp anion 
~ 1 arcmin to the North (|Stavelev-Smith et alj|l9901 ). The 
color profiles are consistent with a mixed stellar population 
contaminated by nebular emission but with no detectable 
dust effects. Due to the inside-out distribution of the burst 
regions the color-color diagrams show increasing nebular 
emission contamination at larger radii. The irregular mor- 
phology further prevents us from separating the burst from 
the older population, and we cannot confirm the presence 
of three distinct popul ations in this galaxy, as suggested 
bv lWestera et all l|2004h . Our population separation method 
rests on simplified assumptions which naturally do not hold 
for galaxies with peculiar morphology. We do observe a pop- 
ulation of < 50 Myr in the burst region and ~ 300 Myr in 
the central red region, however the latter certainly suffers 
from contamination from the burst, and we cannot sample 
a purely old population dominated region at any radius. We 
classify this galaxy as U,M BCG. 

UM238 

Previous observations have compared this galaxy to UM 133 
because its West starburst knot was th e only one visible, 
givin g the galaxy a cometary appearance l|Telles fc Terlevichl 
1997). Our contour and RGB plots indicate at least six such 
regions, which must be super star cluster complexes. The 
brightest one is to the West, as already observed by the 
previous authors, one is to the North, and the remaining 
four to the South-East an d Ea st of the geometric ce nter. 
Both IWestera etail <|2004h and ICuisinier et~al] l|2006h find 
three distinct populations in this galaxy, young (~ 6.8 Myr), 
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intermediate (~ 10 Myr) and old (> 1 Gyr), with the old 
one clearly dominating the total mass of the galaxy. Due to 
the morphology of the galaxy we detect a population > 1 
Gyr in our color-color diagrams but are unable to impose 
any limits on the age of the star forming regions. There is 
an excess of ~ 0.2 mag in all H — K diagrams, which could 
be due to nebular continuum emission, however, we see only 
very mild nebular emission contribution in B — V vs. V — I 
or V — K diagrams. We classify this galaxy as il BCG. 

UM417 

This is an isolated il, C BCG in the sense that there 
are n o detected companions within 1 Mpc projec t ed dis - 
tance l|Campos-Aguilar fc Moleslll99ll '). [Salzer et all l |l989bh 
classify it as a SS object. The color-color diagrams show a 
large contribution from nebular emission, though the effect 
is negligible outside of the starburst region. These diagrams 
also suggest that the burst is of intermediate (Z ~ 0.008) 
or low (Z ~ 0.004) metallicity, and no older than a few 
Myr, which makes it the youngest burst in the sample. 
We see no clear indication of dust in the same diagrams, 
which is consistent wit h the very low extin ction values 
found in this galaxy by iTerlevich et al.l (|l99ll ), who mea- 
sure H a /Hf) = 3.71. Further, the instantaneous burst as- 
sumption does not work well for this galaxy since the inte- 
grated colors for some regions are too red in B — V for a 
given V — I or V — K, which must indicate prolonged star 
formation since it cannot be due to dust. 



1 — 




1— 1 


1 1 1 


V. 


J 


— 1— - 

n 


|J 1 


n n 



-20 -16 -IS -14 -12 

absolute B luminosity [mag] 



Figure 7. Absolute B magnitude distribution of the sample for 
the composite galaxy (upper panel), only the burst (black con- 
tour, lower panel), and the host galaxy (gray filled spikes, lower 
panel). 



from the behavior of the V — R and V — I profiles. The 
first population dominates the range r = 0-10 arcsec, the 
second r — 10-25 arcsec, and the third beyond r > 25 
ar csec. This is c o nsiste nt with the results for this galaxy 
of IWestera et all (|2004r ), who obtain a best fit of the SED 
with three distinct populations. 



UM619 



UM448 

The star formation in this Wolf-Rayet galaxy is occur- 
ring only in the spiral-like structure visible in the RGB 
plot since those are the only regio n s with detect a ble H a 
emission l|Telles fc Terlevichi lil)97r ). ISalzer et all l|l989bT ) 
classi fy it as a H I IH g alaxy due to its spectral proper- 
ties. iNoeske et al.l l|2003T l classify it as il based on NIR 
images. In our contour plot it is evident that beyond 
(iB ~ 25.5 mag arcsec -2 the isophotes become regular 
and elliptical, albeit somewhat boxy. This would make 
an iE BCG classification more suitable. However, further 
examination of this faint structure shows that we can trace 
it out to an additional ~ 35 arcsec beyond the \ib ~ 25.5 
isophote. Figure IB2I in Appendix IB2I demonstrates this, 
clearly showing a broad tidal arm to the North-West. It 
is visible in all optical filters. In the NIR the extended 
component as a whole is much less prominent, however, 
non-negligible signal is still detected from several areas. 
U — B or B — V vs. V — I colors reveal nebular emission 
contamination from the star formation region, and an old 
(> 5 Gyr) metal-poor (Z < 0.001) stellar population in the 
regular elliptical region beyond /j,b ~ 25.5 mag arcsec -2 . 
We are lacking U and H band observati ons so we cannot get 
a feeli ng for the dust content, however. iLisenfeld fc Ferraral 
(1998) obtain M du3t = 1.4 x 10 4 M Q within the optical disk 
of the galaxy, which is a factor of ~ 1000 lower than the 
H I content inside the same region(Afn i = 1.7 x 1O 7 M0). 
The V — K reddening we observe at small radii must 
then be primarily due to gas, rather than extinction. We 
can deduce the existence of three distinct populations 



This is a DANS g alaxy dSalzer et all I"l989b| ) with no de- 
tected companions ( No eske et all 2001 ). It has an H I mass 
of 6.9 x 10 s M Q (|lmpev et alj|200ll h The contribution from 
nebular emission is moderate to negligible in the measured 
colors, judging by how well they are modeled by a pure stel- 
lar population. There is, however, a ~ 0.2 mag excess in all 
H — K diagrams coming from areas brighter than fis ~ 24 
mag arcsec -2 , likely due to nebular emission. The color- 
color diagrams consistently reveal an old population of > 8 
Gyr. In fact, we do not measure integrated colors indicative 
of a population younger than 1 Gyr for any color combi- 
natio n in any of our radial bins. The SFR is reasonably 
high |Burkholder et al.ll200ll . SFR = O.12M /yr) but the 
star forming regions tracing the spiral arms are not domi- 
nant over the underlying host when the radial bins of inte- 
gration are significantly larger than the angular diameter of 
the star-forming knots in the RGB image. We classify it as 
iE BCG due to its regular outer isophotes. 



12 RESULTS AND DISCUSSION 

Before we delve into the discussion we should mention that 
the total apparent magnitudes in Table[4]are consistent with 
tabulated photometry measurements on NED. The largest 
offsets between our measurements and other authors are 
~ 0.2 mag for IIZw40, MK600, and Tol0341-407. This differ- 
ence is independent of the assumed distance to the targets, 
and is likely due to the small variations in Galactic extinc- 
tion corrections applied by the individual authors, and the 
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Figure 8. Binggeli diagram: Central surface brightness of the 
host vs. the absolute B magnitude of the whole galaxy (circle), 
the burst (star), and the host (diamond). The size of the markers 
reflects the scale length of the host. Markers are color— coded by 
the percentage of total light that is attributed to the burst. The 
approximate locations of dwarfs (the area between the two dot- 
ted lines), spirals (S— shaped marker), and LSB galaxies (hatched 
area) are marked for convenience. 




Figure 9. Upper panel: £?-band host /xo in mag arcsec vs. 
log 10 h r in kpc. For reasons explained in § 112.31 we use the pa- 
rameters obtained from the /jg = 26-28 mag arcsec -2 region. 
Lower panel: Absolute magnitude of the burst vs. absolute mag- 
nitude of the (assumed exponential) host. Marker size in both 
panels reflects the scale length. The colors indicate the relative 
burst contribution to the total (burst+host) luminosity of each 
target. Errorbars are overplotted on each marker. 



fact that the used B filter is not identical, and can on aver- 
age result in differences of ~ 0.15 in integrated magnitude. 
In the case of MK600 the three bright nearby stars may have 
influenced our total luminosity measurement. 



12.1 Notes on isophotal profiles 

The weakness of purely isophotal profiles is that they 
become unreliable at fainter isophotes much sooner than in 
the case of elliptical integration. The way we have chosen 
to obtain the isophotal profiles produces a systematic effect 
clearly visible in all radial surface brightness and color 
profiles involving the reference filter, usually B, or V if 
B is missing. The isophotal profile in the reference filter 
becomes systematically brighter as one reaches deeper into 
the sky, then becomes fainter again as the sky residual peak 
is passed, and eventually starts randomly fluctuating as the 
pixel to pixel sky statistics start dominating at the faintest 
levels. The profiles are truncated before the latter happens, 
as the errors become too large beyond 28 B magnitudes for 
all targets. There is no way, however, to prevent bright sky 
pixels in the range ~ 24.5-28 B magnitudes from entering 
the respective magnitude bins, as there is nothing, other 
than spatial position on the image, to identify them as 
sky instead of galaxy. Naturally, this magnitude range is 
approximate and varies with the depth of the observations, 
but for the vast majority of our images the bulk of the 
sky lies in the interval ~ 24.5-28 B magnitudes, with the 
peak at ~ 26 B magnitudes. To reduce this effect, the 
images used for isophotal profiles were masked out as close 
to the galaxy as possible, without risking removing too 
much of the fainter galaxy isophotes. This, however, has 
another, equally undesirable effect of systematically biasing 
the size of the equivalent radius - beyond a certain surface 
brightness (usually hb > 26 mag arcsec -2 ) each new 
radius artificially decreases as the consecutive isophotes 
run out of pixels, which causes the faint end of the profile 
to decrease steeply in an unphysical fashion. Both of these 
effects are much more pronounced in surface brightness 
or color profiles involving the reference filter because the 
areas defined by each magnitude bin were tailored to those 
frames, and hence the effect of the sky is systematic, 
accumulative and pronounced. The rest of the filters, which 
use isophotal areas defined by the reference, have their sky 
brightness peaks at different wavelengths from the reference 
filter and hence the effect is diluted, albeit still systematic, 
and is completely undetectable in color combinations which 
do not involve the reference filter. 

For the majority of our targets the isophotal profiles trace 
the elliptical integration profiles fairly accurately. If the 
slopes (scale lengths) of these two profiles agree then that is 
a good indication that the star formation has not changed 
the galaxy morphology beyond recognition, i.e. that the host 
galaxy has not been completely transformed by the presence 
of the star forming regions, and it retains its structural prop- 
erties (fio and hr) even after the star formation has been 
switched on. 
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12.2 General properties of the sample 

A dwarf galaxy has a total luminosity fainter than — 18 B 
magnitudes, however, this magnitude cut-off varies from 
author to author and has no physical basis. The general 
consensus is that galaxies fainter than this are guaranteed 
to be dwarfs. With this cut-off approximately half of our 
sample consists of luminous BCGs (Figure [7] upper panel), 
while the other half are blue compact dwarfs (BCDs). We 
refer to both types as BCGs. Separating the "burst" and 
"host" populations in each target, as described in § [HI 
reveals that the total luminosity of galaxies brighter than 
~ — 18 is completely burst dominated and the hosts of all 
targets, except the discordant redshift system Toll457- 
262//, have total luminosities consistent with being dwarfs, 
with M B isk > -18 (Figure lower panel). This is also 
seen in the Binggeli diagram in Figure [8j where the hosts 
of the more massive galaxies are clearly located in the 
indicated dwarf-occupied area, with their corresponding 
burst luminosities constituting more than 80% of the total 
galaxy luminosities. In fact, the burst makes up for the 
dominant percentage of the total light output for most 
targets. Only for HL293B and SBS0335-052West, does 
the host constitute the dominant luminosity contribution. 
These are, however, very compact and faint, so it is unclear 
from the surface brightness profiles in Figure [1] if we 
can sample the old population at all, even at large radii. 
Figure [5J further shows that, with diagonally decreasing host 
scale length from left to right, it is the most extended LSB 
dwarfs that host the most dominant bursts in relative terms 
and the brightest bursts in absolute terms. This can also 
be seen in Figure [9] If all hosts with fj,o > 23 mag arcsec -2 
had dominating relative burst contribution regardless of h r 
then this would simply be due to the fact that the contrast 
between the bright burst and the fainter host is much more 
stark if the latter is a LSB galaxy. As such a relation seems 
to be simultaneously h r - and ^-dependent, however, this 
is not a trivial scaling correlation. The more extended the 
host, the lower the central surface brightness (Figure [51 
upper panel) and while this also usually means the host 
is brighter than the more compact objects, it is not the 
brightest host that has the brightest burst in absolute terms 
and it is not the most extended host that has the highest 
burst contribution in relative terms (Figure (9[ lower panel). 
We explain this in terms of star formation efficiency. When 
they are not undergoing a starburst episode, the more 
extended dwarf hosts spend most of their lives forming 
stars at low efficiencies and hence they retain a larger 
H I reservoir. Once the starburst is induced this larger 
available gas mass contributes to a more vigorous star 
formation. In contrast, compact higher surface brightness 
hosts have higher star formation efficiencies when outside of 
starburst mode and thus would end up with smaller H I gas 
masses, i.e. they are more evolved systems. The lower panel 
of Figure [9j further suggests that there are lines of constant 
relative burst contribution which seem to be predominantly 
defined by the similar scale length of the galaxies, and along 
which both the burst and host luminosities can vary in a 
linear fashion. The sample size is rather small, and we can 
only clearly detect the ~ 60% (yellow-green) and ~ 85% 
(red-orange) lines. There are very few targets with burst 
contribution < 40% and hence no low burst contribution 
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Figure 10. Comparing the structural properties of the host 
galaxies in our BCG sample (orange and black stars) to those of 
dEs, dls, LSB and BCD galaxies from th e literature. dE, d l, BCD 
and some LSB data were obtained from [Papaderos et al ] j200Sh 
and references therein, while the giant LSB galaxies were taken 
from lSpravberrv et al.l 1 119951) and references therein. The orange 
stars indicate structural parameters derived from the physical re- 
gion corresponding to /ig = 24—26 mag arcsec -2 , while the black 
stars are from the fainter fig = 26 — 28 region. Note that h r is 
here in units of pc, not kpc. 



lines, however, this is a selection effect - the galaxies of this 
sample were predominantly chosen to be bright starbursts. 

In Figure \W\ we compare the structural parameters of 
our sample to other types of dwarfs, and also to observa- 
tions of blue compact dwarfs (BCD) fro m other authors. 
The d ata in the figure were taken from iPapaderos et al.l 
l|2008h for the dE, dl, BCD and some LSB galaxies, and 
from Spravberrv et all ij 19951 ) and references therein for 
(giant) LSBs. We have obtained structural parameters from 
two separate outer regions, corresponding to fis =24 — 26 
and fiB = 26 — 28 mag arcsec -2 . For the fiB = 24 — 26 mag 
arcsec -2 region the central surface brightness /j,q for our 
sample is consistent with other BCD measurements from 
the literature, as seen in the upper panel of Figure 1101 If 
we use the fainter fiB = 26 — 28 mag arcsec -2 region we 
instead obtain fj,o consistent with those measured for dE 
and dl galaxies. Our profiles are among the deepest BCG 
profiles available in the current literature. Previous BCG 
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studies, which rarely present data significantly beyond 
the Holmberg radius, claim that for a given total host 
luminosity M^ost, the hosts of BCGs ten d to on average 
have a /io brighter than dEs and dls (e.g. IPapaderos et al.l 
2008). Our data clearly demonstrate that if one samples 
fainter host regions there is no separation between dEs, 
dls and BCGs in terms of /io vs. Mhost- A similar trend is 
observed when we compare the scale length derived from 
(is = 24 — 26 and /ib = 26 — 28 mag arcsec -2 to literature 
data. The brighter region gives h r values consistent with 
the BCD occupied region in the lower panel of Figure 1101 
although the average h r is slightly larger for our sample 
even in this region. The fainter (is = 26 — 28 mag arcsec -2 
region moves the galaxies to scale lengths which are 
inconsistent with the BCD observations from the literature. 
While some data points are now clearly in the region 
occupied by dEs and dls, quite a few have h r even larger 
than that, with \ogh r > 3.5. This is not a problem since 
the dE and dl observations themselves rarely probe surface 
brightness levels as faint as our /xb = 26 — 28 mag arcsec -2 
region. While the dE profiles are unlikely to change slope 
in deeper observations, this is entirely possible for many dls. 

Another possibility exis ts, however, when we co nsider the 
giant LSB galaxies from ISpravberrv et al.l |l995). Both our 
extreme scale lengths (log h r > 3.5) and their corresponding 
central surface brightnesses (fio > 24 mag arcsec -2 ) are 
perfectly consistent with the structural parameters of these 
giant LSB galaxies (LSBG). There is no convention on the 
definition of LSBG, although usually galaxies with central 
surface brightness fainter than /ib ~ 23 mag acrsec -2 
and more luminous than M B ot ~ — 14 mag are considered 
to fall in this category. The location of our galaxies in 
Figure [10] suggests a tentative connection between LSBG 
and ( at least some) BCG hosts. This result is consistent 
with iBergvall et ail <| 19991 ). who argue that luminous 
starbursts in BCGs are fed by LSBGs and not dls due 
to luminosity, metallicity and H I mass considerations. 
Additionally, this result further exemplifies the importance 
of sampling the host well away from the effects of the 
burst. While this may vary in individual cases, the surface 
brightness profile of BCGs seems significantly affected by 
the burst almost down to the Holmberg radius. However, 
while a difference in burst contribution would explain the 
varying /io and h r between the two regions, it cannot be 
excluded that the physical structure is indeed changing at 
the faintest levels, and that many luminous BCG s have 
double-disk profiles (e.g. iHunter fc Elmegreenll2006l ). 

The majority of the BCGs show color gradients in purely 
optical (B — V), purely NIR (H — K), and optical-NIR 
hybrid (V — K) colors, and flat profiles at large radii. In 
terms of the luminosity the highest burst contributions 
are found in galaxies which light distribution we can 
trace down to (ib ~ 28 mag arcsec -2 , and for which we 
observe a clear break in this distribution (i.e. a change 
in stellar population) beyond (ib > 26 mag arcsec -2 . 
These are ES0185-13 (B band burst contribution 95%), 
ES0249-31 (86%), ESO400-43 (94% in V), HE2-10 (89%), 
Tol0341-407 (83%), and UM448 (88%). The burst domi- 
nance is not solely due to defining everything above this 
new very faint host as a burst. A profile with constant 



slope in the range fis ~ 24.0-28.0 would indicate that 
there is no change in population over that region and 
since the radial distances are large, presumably this single 
population is the host. Such targets also display varying 
burst-domination percentages. Some are rather high, e.g. 
SBS335-052East (86%), MK930 and ESO338-04 (both 
with 94%), etc., while others are much less, e.g. UM238 
(60%), or host dominated instead, e.g. HL293B (16%) and 
SBS335-052West (39%). Note that previous studies rarely 
present the surface brightness profile much further than 
the Holmberg radius with comparably sma ll errors (e.g. 
Papaderos et al.l Il996l ; iDoublier et alj 1 19991 ; ICairos et alj 
2001al ). so for targets in our sample which have a profile 
break around [is ~ 26.0 what we here define as host is 
a previously undetected population of different structural 
parameters and photometric properties than what is 
measured at the smaller radii between /is ~ 24-26 mag 
arcsec -2 . For ES0185-13 and UM448 the now revealed host 
is a true LSB galaxy, with fj, = 24.7 ± 0.12, 24.33 ± 0.09 
mag arcsec -2 , and h r = 5.9±0.3, 9.3±0.4 kpc, respectively. 



12.3 Color— color diagrams 

The parameters that can influence the colors of a galaxy 
in the optical and NIR regimes are the metallicity, the age 
of the stellar population, the internal extinction, the star 
formation history, and the nebular component, be it ionized 
in situ or by photons "leaking" from nearby star forming 
regions. These can conspire to either enhance or counteract 
each other's influence on the galaxy colors, which is why the 
best way to analyze the colors we have obtained for various 
physical areas and stellar populations is to use spectral 
evolutionary models (SEMs) based on stellar population 
synthesis. We have chosen to present only B — V vs. V — K 
because this diagram helps break the metallicity degeneracy 
more than the rest. 

In Figure [TT] we compare our data to the two cases of 
a pure stellar population (right panel), and of stellar 
continuum contaminated by nebular emission (left panel). 
The data we plot, in descending order, are the total color 
of each galaxy measured down to the Holmberg radius from 
Table |4j the estimated color of the burst from Table [7] 
and the color of the host regions defined over two radial 
ranges corresponding to /j,b = 24-26 and /j,b = 26-28 mag 
arcsec -2 from Table \5\ (see § [7] for details). The errorbars 
in the individual subplots in this figure are the same as in 
the corresponding tables. Note that for the burst colors all 
errorbars are underestimated, as they only contain the very 
small numerical error of integrating above the interpolated 
exponential disk, and the zero point uncertainties of the 
respective nights. 

The models with nebular emission are from the Yddrasil 
spectral synthesis code, as already described in § 1111 While 
these are suitable for the young population and for regions 
with nebular emission contribution, the analysis of NIR 
fluxes and colors is very sensitive to the treatment of 
thermally pulsing AGB stars. We therefore also present 
an ind ependent set of models based on the iMarigo et al.l 
(2008) isochrones, which feature one of the most up-to-date 
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Figure 11. Integrated optical/NIR colors with SEMs tracks with (left panel) and without (right panel) nebular emission. From top to 
bottom in both panels the plotted data points are the total colors over the entire area of the galaxy, the estimated burst colors, the colors 
of the /ig ~ 24 — 26 mag arcsec -2 region, and the colors of the fig ~ 26 — 28 mag arcsec -2 region. The data are corrected for Galactic 
extinction. The errorbars include different uncertainties depending on the case (see § 112,31 for details). The number of targets in the plots 
is not equal since the burst/host separation fails completely for a number of objects with irregular morphology. The tracks in the left 
panel include nebular emission and assume an instantaneous burst, with metallicities Z = 0.0004 (dotted), Z = 0.004 (dash-dotted), 
Z = 0.008 (dashed), and Z = 0.02 (line). The ages of 10 Myr (open square), 1 Gyr (open circle), and 14 Gyrs (open triangle) are marked 
for convenience. In the right panel the tracks are for a pure stellar population with a Salpeter IMF, M m j n = O.OSMq, M max = 120Mq, 
an e-folding time of 10 9 yr, and with metallicity Z = 0.001 (dotted), Z = 0.004 (dash-dotted), and Z = 0.02 (line). The ages of 1 (open 
circle), 3 (open left triangle), 5 (open square), and 14 Gyrs (open triangle) are marked for convenience. 
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treatments of this evolutionary stage currently available. 
However, as these models lack nebular emission, they are 
only suitable for the outskirts of our target galaxies, where 
nebular emission tends to be subdominant. These model 
tracks assume a Salpeter initial mass function (IMF) and a 
star formation decay rate of 1 Gyr, suitable for early-type 
galaxies. The star formation history (SFH) is introduced 
by weighting each consecutive point with the average star 
formation of the previous points, and is suitable for an old 
stellar population with no recent burst and very low SFR. 
The predictions of both SEMs tracks are shown for several 
metallicities ranging from very metal-poor [Z — 0.001 or 
Z = 0.0004) to solar (Z = 0.02). 



individual galaxies the contamination may be negligible 
depending on the morphology, location and extent of the 
burst regions, and the distribution of the gas over the disk 
of the galaxy, but this would have to be examined on a 
case to case basis before any measurements over this region 
are considered reliable. This is not possible if the data are 
not deep enough to reach beyond the Holmberg radius with 
small errors. Our data allows this for the majority of the 
targets. We have therefore adopted the structural param- 
eters of the host derived from the clearly host-dominated 
(iB — 26-28 mag arcsec -2 region when analyzing possible 
correlations between the structural parameters and other 
galaxy properties. 



All SEMs tracks plotted in Figure [TT1 assume zero redshift, 
which is a coarse approximation for our case since half of 
our sample is at a redshift of ~ 0.02. The B — V colors of 
the young component can vary by as much as 0.3 mag from 
redshift 0.00 to redshift 0.02. Further, even though we have 
used "standard" BVK filters in each observing run, filters 
over similar wavelength ranges from different observatories 
show small differences in transmission and shape. These 
differences can amount to e.g. ~ 0.15 mag for the B band. 
Summa summarum, the tracks plotted in Figure QTJ are 
only intended to guide the eye and give an overall estimate 
of the age and metallicity of the data. A much more 
careful and detailed analysis of all colors, with SEMs tracks 
tailored to each galaxy and filter, lies beyond the scop e 
of this paper and will be provided in lMicheva et al.1 (|2012bh . 

In the B — V vs. V — K plane, contribution from nebular 
emission to the stellar continuum shifts the predicted colors 
redwards in B — V and bluewards in V — K, an effect 
which is much more drastic for ages younger than 1 Gyr 
(left panel of Figure ITT)) . In fact, ignoring nebular emission 
(right panel of Figure 1 1 1 p leads to clearly irreconcilable 
differences between the observed and the predicted colors 
for the galaxy as a whole, for the burst, and also for 
the (lb = 24-26 mag arcsec -2 region. Even though we 
present both sets of tracks, with and without nebular 
emission, it is clear that e.g. the burst colors can only 
be explained by the Yggdrasil tracks because the young 
regions of the galaxies have larger amount of ionized gas 
and more young stars present to ionize the gas. In contrast, 
for the fiB = 26-28 mag arcsec -2 region which contains 
predominantly old stars and much less gas, the colors are 
best reproduced by the stellar population tracks seen in 
the right panel. Not all targets have detectable K band 
flux at these radial distances even given the large bin size 
of the integration, and the ones that do have quite large 
errorbars, but the trend towards an old stellar population 
is nevertheless clear. We note that the (is = 24-26 mag 
arcsec -2 region, which is at moderately large radii from the 
burst, still shows heavy nebular emission contamination. 
The colors here are on average similar to the total galaxy 
colors and much bluer in V — K than in the (is ~ 26-28 
mag arcsec -2 region. This indicates that there is still 
significant contribution from the young population to the 
light distribution of the galaxy. This region is therefore 
not to be relied upon for the measurement of structural 
parameters of the host such as the scale length ha and 
the central surface brightness (iq. Undoubtedly, for some 



In the \ib = 26-28 mag arcsec -2 region there are a number 
of galaxies with measured V — K which appear to be too 
red and indicate super-solar metallicities. The errorbars 
are large and all except the two most extreme points can 
be reconciled with the pure stellar population tracks with 
lower than solar metallicities, but we will nevertheless 
examine these on a case to case basis. Solar or near-solar 
metallicities are unusual for such large radii since we expect 
the (stellar) metal l icity to drop with increased radius from 
the center (lBekki|[2008l ). The latter only holds, however, 
for regular morphology galaxies. The deviating galaxies, 
in order of decreasing V — K excess, are UM619, UM238, 
UM160, and ESO421-02, which all cluster at the solar 
metallicity track. Statistically one would expect outliers like 
this in photometry measurements without this having to 
mean that there is something wrong with the measurements 
themselves. U M619 has nearly solar m etallicity according to 
the literature l|Burkholder et al"1l200ll ). and most of its star 
formation is distributed along the spiral arms at significant 
distances from the center. UM238 has all of its bright star- 
forming cluster complexes in the outskirts of the host, some 
of which fall inside the integration bin corresponding to the 
(iB = 26-28 mag arcsec -2 region. Similarly, UM160 has 
highly irregular morphology, with many blue regions found 
inside of the relevant integration radial range. ESO421-02, 
on the other hand, has outskirts clearly dominated by a red 
host, but this galaxy shows nearly constant V — K color for 
all physical regions we examined - the difference between 
the total galaxy color and the color of the two outer regions 
is A(V — K) < 0.06 mag. Thus all of these measurements 
are either due to the outer regions being contaminated by 
starburst contribution, increasing its apparent metallicity, 
or the galaxy as a whole shows no metallicity gradient. 
The measurements contaminated by burst are therefore 
better compared to the Yggdrasil tracks, which include 
nebular emission. Specifically the Z — 0.008 track with 
an age < 1 Gyr seems suitable for most of these data points. 

Previous studies of BCGs have detected host V — K col- 
ors which are difficult to reconc ile with a norma l stel- 
lar population of low metallicity feergvall fc Ostlinl |2002| ; 
iBergvall et"ai1l2005l ; IZackrisson et al.1 12006). Some extreme 
cases, such as Haro 11, were impossible to reconcile with 
any metallicity. However, judging by the B — V vs. V — K 
diagram, the case for an abnormally red host (a.k.a. the 
"red halo" phenomenon) in BCGs seems rather weak for this 
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Table 8. Minimum Potrosian Ap asymmetries measured in each filter over the area enclosed by the Petrosian radius 
r[r/(0.2)], here given in kpc for each galaxy and filter. 
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sample. Though there are only eight color measurements in 
the hb = 26-28 mag arcsec -2 region, and the errorbars are 
large, we see no clustering of the data points around high 
metallicity tracks or a general tendency towards abnormally 
red V — K colors. We have also verified this with a B — V vs. 
V — I diag ram (not p r esent ed here) using the same two mod- 
els. In lMicheva et al.l l|20T0t ) we further ref uted the extremely 
red color claimed for the h ost of Haro 11 l|Bergvall fc Ostlm] 
120021 ; iBergvall et al.ll2005l . V - K ~ 4.2 ± 0.8), and mea- 
sured instead the more modest color of V — K ~ 2.3±0.2 
over the same physical region. In other words, the present 
work casts doubt on whether any BCG truly has an ab- 
normally red host, since Haro 11, previously the most ex- 
treme case of a red BCG host, is no longer extreme, and 
we find no clear indication of abnormally red colors in the 
present sample. Though there are currently no detections of 
an abnormally red host in individual BCGs left undisputed, 
here we cannot decisively conclude that the hosts of BCGs 
have normal co lors. We will make suc h a determination in 
a future paper iMicheva et al. I (|2012bh . Note that the lack 
of abnormal colors in BCG hosts has no direct bearing on 
the extremely red colors of halos aroun d high- and low sur- 
face brightness stac ked edge-on disks (|Zibetti et alj |2004| ; 
IBergvall et aill2010l ). 

While we have noted that the fiB = 24 — 26 mag arcsec -2 
region for this sample contains a strong contribution from 
the star forming regions, from Figure [TTJ we can still make 
the observation that all host galaxies in the sample do not 
appear to be younger than 5 Gyr. The same conclusion, 
albeit with fewer available data points and larger errorbars, 
can be drawn from the fiB = 26 — 28 mag arcsec -2 region 
in the same figure. The implied redshift at which the bulk 
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Figure 12. Histograms of Ap Petrosian (left) and A' H Holmberg 
(right, black) asymmetries. Overplotted in the right panel (red) 
is the dynamical Ad yn asymmetry. Bin size of 0.05 is used. The 
plotted data are from Tables [8] and [9] 
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of the star formation in the host must have occurred is thus 
z ~ 0.5 (8 to 9 Gyr after the Big Bang). While older stars 
could of course be present in all of our targets and can 
easily remain unseen in these data, it appears unlikely that 
the prototypical host for this sample could have been a red 
and dead elliptical galaxy since those have colors consistent 
with the Hubble age (compare to KIG732 in Appendix [X] 
with B -V ~1, V - if ~ 3 at large radii). 



12.4 Asymmetry and concentration 

The histograms in Figure fTS] show the distribution of Pet- 
rosian asymmetries (left panel), and "Holmberg" asymme- 
tries (right panel) of the galaxies in the sample. In the case of 
the Petrosian asymmetry the size of the area over which it is 
measured is small (see Table|8}, and the faintest isophote in- 
cluded in that area is often as bright as fi ~ 21 mag arcsec -2 . 
The Petrosian asymmetry is therefore mainly flocculent in 
nature - it is a measure of the deviation from the symmet- 
ric ground state due to the presence of bright SF regions. 
As an example of how flocculent asymmetry can increase 
the total asymmetry of an otherwise symmetric galaxy, con- 
sider the case of HL293B. This galaxy has regular elliptical 
outer isophotes and should thus merit a rather small asym- 
metry, but instead has Ap ~ 0.3 due to its off-center star 
formation region. One would therefore expect the composite 
asymmetry to correlate with relative burst contribution in 
the sense that as the burst contribution decreases the dom- 
inating asymmetry component changes from flocculent to 
dynamical, where the latter is a measure of low frequency 
structure. This does indeed seem to be the case for about 
half the sample as seen in Figure 1131 As the contrast be- 
tween star forming regions and faint dynamical features de- 
creases with the burst contribution, the measured compos- 
ite asymmetry increases as more weight is placed on the 
dynamical morphology. There appears to be an upper limit 
to the burst contribution at about > 90%, above which the 
sampled asymmetry is purely flocculent. Perhaps not sur- 
prisingly, the galaxies that fall within this category are the 
BCGs consistent with true LSB hosts, as identified by Fig- 
ure [TO] and discussed in § 112.21 The contrast between the 
bright burst and the LSB host is here very high, which 
would make the effect of dynamical asymmetry negligible 
since many tidal features or extended irregular structures 
are too faint to fall inside of the Petrosian r[7j(0.2)] radius. 

The dynamical asymmetry should here be better repre- 
sented by what we call the Holmberg asymmetry (Table [9j. 
The term is a bit misleading since the Holmberg radius is 
only defined in blue passbands, and we use it purely for ab- 
breviation purposes. To reduce the effect of the bright, often 
irregular SF regions on the total asymmetry budget, and to 
compensate for the difference in distances to the galaxies 
in the sample, the images have been boxcar smoothed by a 
filter of size lxl kpc. This has the effect that predominantly 
large scale structures remain and the measured asymmetry 
should therefore be mainly dynamical, not flocculent. Our 
hope that the asymmetry obtained in such a fashion would 
be a reasonable proxy for the dynamical asymmetry was 
crushed by our findings in § 112.21 and § 112.31 where we 
established that for this sample of luminous BCGs there 
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Figure 13. Relative burst contribution to the total luminosity 
in percent vs. the B band Petrosian Ap asymmetry for the sam- 
ple. The filled black circles are the galaxies with true LSB hosts, 
flo ^ 23 mag arcsec -2 and h r > 3 kpc. The open circles are the 
remaining galaxies in the sample. The size of the markers reflects 
the respective scale length. 
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Figure 14. Petrosian asymmetry (Ap) vs. all alternative asym- 
metry measurements. Ap is the Holmberg asymmetry without 
any smoothing of the images, while a 1 X 1 kpc smoothing box 
has been applied to A' p (Petrosian), A' H (Holmberg), and A^ yn 
(dynamical). The uninterrupted lines are least square fits to the 
similarly colored data points. The dashed line has a slop e of 1. The 
dotted line is an extrapolation from the correlation of IConselicei 
||2003|) for E, SO, Sa-b, Sc-d and Irr, which lie beyond A P ~ 0.6. 
All asymmetries here are f r om th e I band, in order to facilitate 
comparison with lConseilcel <2003h . 



is still a strong contribution from the starburst down to 
the \ib ~ 26 mag arcsec -2 isophotal level, which is nearly 
at the Holmberg radius. We will denote the Petrosian 
asymmetry with Ap, the unsmoothed Holmberg asymmetry 
with Ah, and the smoothed Holmberg asymmetry with 
A' H , where the prime indicates smoothing. In Figure [14] 
we find a positive correlation between the Petrosian Ap 
and Holmberg Ah, A' h asymmetries in the sense that for 
a given galaxy large/small Petrosian nearly always implies 
large/small Holmberg asymmetries. This makes sense for 
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starbursting luminous BCGs, as one would expect there 
to be a strong positive correlation between flocculent and 
dynamical asymmetries. A (wet) merger would lead to star 
formation, and the former would increase the dynamical 
component, while the latter would increase the flocculent 
component in the composite asymmetry. In our case, 
however, there are indications that the Petrosian and 
Holmberg asymmetries are measuring the same (flocculent) 
asymmetry component in over half the sample, with the 
numerically different Holmberg values being solely due 
to the applied smoothing, and not to the sampling of a 
dynamically dominated asymmetry. For a few cases (open 
circles in Figure II3[) the Holmberg asymmetry successfully 
measures predominantly the dynamical component, as 
exemplified by HL293B. From a Petrosian asymmetry 
of Ap(B) ~ 0.3 the value has now dropped to ~ 0.06, 
which is much more suitable for a galaxy with such regular 
morphology. 

To obtain a more universally successful measure of the 
dynamical asymmetry for our sample, we tweaked the 
Holmberg A' H asymmetry to allow it to ignore star forming 
areas. All regions brighter than /i = 25 mag arcsec -2 in the 
optical, and fi = 21 mag arcsec -2 in the NIR were set to a 
constant value of 25 and 21 mag arcsec -2 respectively in 
the two wavelength regimes. This effectively circumvents 
the contribution of the flocculent asymmetry component 
to the total asymmetry budget, and it has the effect that 
more weight is given to the fainter structures. This modified 
asymmetry is a better estimate of the dynamical asymmetry 
component, and is given in Table [9] We will refer to this 
asymmetry as simply dynamical asymmetry, Adyn- For 
targets for which the Holmberg asymmetry (A' H ) already 
measures the dynamical component, the Adyn values are 
very similar (e.g. HL293B, MK996, and other ellipticals). 
For targets for which the flocculent asymmetry dominated 
even the Holmberg values, Adyn varies from case to case 
depending on the morphology. ES0185-13, for example, 
now shows its more asymmetric nature, as the decreased 
flux contrast gives its faint tidal tail a chance to contribute 
significantly to the total asymmetry budget in the optical, 
resulting in Ad yn (B) ~ 0.2 instead of A' H (B) ~ 0.1. In the 
NIR the values stay the same, A(NIR) ~ 0.1, since the 
tidal tail is invisible at these wavelengths. 

Figure [14] summarizes the behavior of all of these al- 
ternative asymmetry measurements versus the standard 
Petrosian asymmetry. In order to illustrate the effect of 
the resolution, both a smoothed (A') and unsmoothed 
(A) version of the Petrosian (Ap, A' P ) and Holmberg 
(Ah, A' h ) asymmetries are included in the figure. The 
correlation between Petrosian and Holmberg asymmetries 
(both smoothed and unsmoothed) is clear, with Pearson's 
R — 0.9; 0.8; 0.8, and somewhat similar slopes for the 
fitted lines m = 0.66; .50; .49 for A P vs. A H , A' H , A' P , 
respectively. IConselicd (|2003l ) finds m = 0.67 for spheroids 
(E,S0), early and late type disks (Sa-b,Sc-d), and irregulars 
(Irr). From all of our asymmetry measurements, A' H is 
best suited for direct comparison with Concelice's "global" 
A G (R) asymmetry. We get A' H {I) = 0.50 x A P (I) + 0.03. 
This is consistent with Figure [13] where for ~ 50% of the 
targets the total asymmetry is completely dominated by 



the star formation. 

In contrast to our previous comparison, in Figure [14] we do 
not see any correlation between the Petrosian (flocculent, 
Ap) and the dynamical Adyn asymmetries. This makes 
sense because by substituting all star forming regions with 
a constant brightness we are effectively only sampling 
the dynamical asymmetry component of the host, and 
not of the composite galaxy. One does not expect any 
particular dependence between the flocculent asymmetry 
of the composite galaxy and the dynamical asymmetry of 
the host, except in some merger cases. We say some merger 
cases, because there are also examples of very symmetric 
mergers, such as ESO421-02. 

For a given galaxy, the Petrosian asymmetries do not 
show much variation from filter to filter, and preserve 
similar values even when going from the optical to the 
NIR regime (Table [5}. The same applies to the Holmberg 
asymmetries, which show similar negligible change with 
filter for the majority of the galaxies (Table [9j. The NIR 
asymmetry, regardless of how it is measured, should be 
more sensitive to the morphology of the old stellar popula- 
tion, and hence to the dynamical contribution to the total 
asymmetry. On the other hand, the optical asymmetry 
is more attuned to the flocculent asymmetry component 
since the young stellar population and the nebular gas 
would dominate the optical emission. The near constancy 
of the asymmetry across filters can then only imply one of 
two things - that both the optical and NIR regimes are 
probing predominantly the flocculent asymmetry, which 
here dominates even in the NIR; alternatively, that the 
flocculent and dynamical asymmetries are equal and there 
is a perfect anti-correlation between the filter sensitivity 
towards the two asymmetries, in the sense that e.g. the B 
band is very sensitive to the flocculent asymmetry and very 
insensitive to the dynamical asymmetry, and that this trend 
is reversed as one moves to the NIR. The latter is unlikely 
since it would require a remarkable degree of fine-tuning. 
Additionally, even if a specific filter is more sensitive to 
a particular population, it is not completely insensitive 
to the other population. Therefore when it comes to the 
asymmetry of the composite galaxy, i.e. either Petrosian 
or Holmberg, then both the optical and NIR asymmetries 
must be flocculent in nature. Note, however, that while 
the dominance of the flocculent asymmetry in both the 
optical and the NIR may seem intuitive and self-evident, we 
find large variations in asymmetry values between the two 
wavelength regimes for the volume-limited sample of emis- 
sion line galaxies in I Micheva et al] l|2012al ). which correlate 
with morphological class. The constancy of the Petrosian 
and Holmberg asymmetry for this sample is due to the 
selection criterion which favors very luminous BCGs with 
intense star formation. These are predominantly burst dom- 
inated, as seen from Table [7] and Figure PHA with more than 
half the sample showing a relative burst contribution > 70%. 

In Figure [15] we compare the concentration from Table [10] 
to the Petrosian Ap asymmetry for our sample. Included 
in the plot are the spheroids (E and SO), early (Sa-b) 
and late type (Sc-d) dis ks, and irregular galaxies (Irr) 
from lConselice et all \20fXt) and references therein. The nor- 
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Figure 15. Concentratio n vs Petrosian Ap As ymmetry for our sample (red text) compared to spheroids, early and late type disks 
and irregulars taken from IConselice et al.l teOOOh . Abbreviations stand for: £JS5=ES0185-13; £2,££>=ESC>249-31; £33S=ESO338-04; 
E400=ESO4OO-43; E421 =ESO421-02; £^62=ESO462-20; H210=HE2-10; HL=HL293B; Zw=UZw40; M600=MK6OO; M000=MK9OO 
M030=MK93O; M000=MK996; S335S=SBS0335-052E; 5335VF=SBS0335-052 W; T341 =TOL0341-407; T1457II=TOL1457-262II 
T^57/=TOL1457-262I; «3=UM133; J60=UM16O; 23S=UM238; ^7=UM417; 4^S=UM448; «0=UM619. 



mal galaxies exhibit a concentration-asymmetry correlation 
in the sense that the very symmetric galaxies are highly con- 
centrated, with spheroids topping the scales in this correla- 
tion. Our sample of luminous BCGs occupies the opposite 
end of the extreme - on average they display higher asym- 
metries than the other galaxy types, and are simultaneously 
much less concentrated. Although we have plotted B band 
concentration and asymmetry, this trend is qualitatively pre- 
served in the rest of the filters (Tables 151 and [TIT)) . The only 
outlier in this relation is SBS0335-052E. As already dis- 
cussed in Ulll this galaxy has an extremely compact central 
region with fairly regular and bright outer isophotes, which 
would bring down its flocculent asymmetry, but is compara- 
tively very extended in the optical, which would bring down 
its concentration. 



13 CONCLUSIONS 

We have obtained and analyzed very deep UBVRIHKs 
broadband imaging data for a sample of 24 blue compact 
galaxies (BCGs). These are a part of a larger data set 
consisting of 46 BCGs, which comprises the most extensive 
study to date of extremely deep optical and NIR surface 
photometry of this type of galaxies. The remaining galaxies 
are a volum e limited sample and are therefore analyzed 
separately in lMicheva et al] |2012al ). 

In this work we have presented isophotal and elliptical in- 
tegration surface brightness and color profiles, as well as 
extremely deep contour maps and RGB images. We have 
measured the total integrated colors of the targets, as well 
as over different regions of the host population. Where the 



morphology allows, we have estimated the burst luminos- 
ity and colors. All measurements were compared to stellar 
evolutionary models both with and without contribution by 
nebular emission. Structural parameters of the host popula- 
tion were derived from a least-squares fit of an exponential 
disk to the outskirts of the light distribution profiles. Our 
main results can be summarized as follows. 

• For the targets lacking a morphological classi fication we 
have p rovided such a classification following lLoose fc Thuanl 
l|l986l ) based on an analysis of our contour and RGB plots. 
Such a classification is highly dependent on which faint sur- 
face brightness level one is looking at, since the structure in 
the outskirts can significantly change its physical properties 
beyond the Holmberg radius. 

• We find no significant difference in the shape and slope 
of isophotal and elliptical integration surface brightness and 
color profiles in the overlap region, except in the cases where 
the burst is not centrally located. 

• Our deep data reveal a previously undetected low sur- 
face brightness component beyond /ib = 26 mag arcsec -2 
for ES0185-13 and UM448. This new low surface brightness 
component has fio = 24.7±0.12, 24.33±0.09 mag arcsec -2 , 
and h r — 5.9 ± 0.3, 9.3 ± 0.4 kpc, respectively. We note a 
similar low surface brightness structure for UM133 beyond 
Hb = 28 mag arcsec -2 . Even though this is at the edge of 
our reliability limit for the B band, the deep contour images 
in the Appendix indicate that this feature could very well 
be real, but we have not dwelled on it further. 

• Regardless of the total luminosity of the BCG, the un- 
derlying host is consistent with being a dwarf with Mb ^ 
— 18, and in the majority of the cases it also qualifies as a 
true low surface brightness galaxy (LSBG) with no 23 
mag arcsec -2 . 
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Table 9. Minimum asymmetries measured in each filter. The 
two numbers per filter per galaxy are Holmberg A' H asymme- 
tries measured over the area enclosed by the Holmberg radius 
r(fj, = 26.5) in the optical and by r(/i = 23) in the NIR (top 
value), and the dynamical A^yn asymmetries, with regions 
fj, < 25 (p, < 21) set to 25 (21) mag arcsec" 2 in the optical 
(NIR) (bottom value). The images are pre— processed by a box- 
car average of size lxl kpc. Dots indicate that the measured 
value was unreliable due to severe masking effects. 
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• The most extended low surface brightness hosts have 
the brightest bursts in absolute terms, and the most domi- 
nant bursts in terms of percentage of total B band luminos- 
ity. 

• We detect an optical bridge between the companion 
galaxies ESO400-43A&B which we can trace down to the 
[iv ~ 28th mag arcsec -2 isophotal level. In contrast, we de- 
tect no optical bridge between SBS0335-052E&W at surface 
brightness levels brighter than /is ~ 28 mag arcsec -2 . 

• The outer regions of a BCG do not always represent a 
region that is host -dominated, even if they have vaguely el- 
liptical isophotes. Derivation of structural parameters there- 
fore either needs to be performed on a case to case basis, 
which will complicate the analysis due to the inherent in- 
consistency of such a strategy, or the interpretation needs 
to allow for the obviously different structural properties of 
BCGs. Specifically, even in the cases of centrally-located 
star-forming regions we find the region defined by \ib = 24- 
26 mag arcsec -2 to be on average strongly contaminated by 
nebular emission, which will affect the slope of the light dis- 
tribution and hence both the derived scale length and the 
central surface brightness attributed to the host. 

• The fainter fiB = 26-28 mag arcsec -2 region is on aver- 



Table 10. Concentration parameter for each filter as 
defined by Equation [2] 
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age free of nebular emission. Using this region for the deriva- 
tion of the structural parameters, we find a significant differ- 
ence between fio and h r for this region and the \xb — 24-26 
mag arcsec -2 region. With the parameters obtained from 
the fainter region we find no difference between the hosts of 
our luminous BCGs and other types of dwarfs like dEs and 
dls. For the truly starbursting BCGs in the sample we find 
evidence suggestive of a tentative connection between low 
surface brightness galaxies and the BCG hosts, since they 
have comparable central surface brightness, scale lengths, 
and total luminosities. 

• We have achieved successful separation of the floccu- 
lent and dynamical asymmetry components for our galaxies 
using a very simple idea. There is no correlation between 
the dynamical and the flocculent asymmetry components 
(Pearson's R = 0.03). 

• In contrast to spheroids (E,S0) and disks (Sa-b, Sc-d), 
luminous BCGs with LSBGs as hosts have total asymmetries 
completely dominated by the star formation event, i.e. by 
the flocculent component. This trend persists in both the 
optical and NIR regimes. 

• In the concentration vs. asymmetry parameter space 
the BCGs are clearly distinguished from spheroids (E,S0) 
and disk galaxies (Sa-b, Sc-d), occupying a region with high 
asymmetries and low concentrations, as expected. 
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APPENDIX A: KIG 732 

This is a textbook pure elliptical galaxy. It is not a part of 
the sample, but it was nonetheless observed by us, and we 
present it here since the observations are somewhat deep and 
previously unpublished. The exposure times are 20 minutes 
each for BVI, and 25 minutes for K and were obtained dur- 
ing 2005 at the NOT. Note that there were only two I band 
images taken during the corresponding observing night. We 
were therefore unable to perform a fringe correction. The 
effect of the lingering fringes is visible in the V — I radial 
color profile beyond r ~ 40 arcsec. 

APPENDIX B: EXTREMELY FAINT 
CONTOURS 

Bl ES0185-13 

In Figure [Bll we investigate the faint features of this galaxy. 
The faint structures are undetected in the NIR. 

B2 UM448 

In Figure lB2l we investigate the faint features of this galaxy 
in all filters. Note that the faint extended features are visible 
even in the K band. 

B3 Extreme contours 

In Figure [B3l we show extremely deep contours for a selec- 
tion of galaxies from our sample. The galaxies not shown 
either reveal no new information compared to the shal- 
lower contours in Figure Q] or deeper contours could not be 
reached. The images are all B band, except for ESO400- 
43A&B where we used the V band. In units of mag arcsec -2 
the contour levels are: 

• ES0185-13: red - 25.8, black - 27.8 

• ESO338-04: red - 25.7, black - 27.7. Note the extended 
structure to the South- West. 

• ESO400-43: red - 25.9, black - 27.9. Note that there is 
an optical bridge between the two companions. ESO400-43B 
is not included in our sample since it was constantly at the 
edge of the frames. 



• ESO421-02: red - 25.9, black - 27.9. The regular shape 
of the outer contour is further indication that the two dwarf 
components involved in this merger could not have under- 
gone a fly-through since there is no detectable tidal debris. 

• HE2-10: red - 25.5, black - 27.2 

• HL293B: red - 25.7, black - 27.7 

• IIZw40: red - 23.5, black - 26.0 

• MK600: red - 25.8, black - 27.8. The 25.8th contour 
is the last one we can close around the galaxy. At fainter 
isophotes the contours merge with those from the bright 
surrounding stars. 

• MK900: red - 25.8, black - 27.8 

• MK930: red - 25.5, black - 27.5. The features to the 
North-East are possibly tidal debris from the recent merger. 

• MK996: red - 25.8, black - 27.8 

• SBS0335-052: red - 25.8, black - 27.8. Note the absence 
of an optical bridge between the East and West companions. 

• Tol0341-407: red - 26.0, black - 28.0. At faint isophotes 
the galaxy radically changes shape, looking more elliptical 
and regular. 

• UM133: red - 27.9, black - 29.1. These are the deepest 
contours for the sample and they are consistent with the 
feature we see in the surface brightness profile beyond [is ~ 
28 mag arcsec -2 . The bright star to the South- West of the 
galaxy goes below B = 29 mag within 40 arcsec, so it is not 
contributing to the structures visible to the East beyond the 
red contour. 

• UM160: red - 26.0, black - 28.0. The galaxy seemingly 
completely disappears beyond the 28th isophote. Since it 
appears to be isolated, this could be a projection effect. 

• UM238: red - 26.0, black - 28.0. The structure to the 
South- West belongs to a bright star. It does not interfere 
with any measurements for this galaxy inwards of the 28th 
isophote. 

• UM417: red - 26.0, black - 28.0 

• UM619: red - 26.0, black - 28.0 
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Figure Al. KIG732. Left panel: Surface brightness and color radial profiles for elliptical (open circles) and 
isophotal (red circles) integration. Upper right panel: contour plot based on the B band. Isophotes fainter than 22.0, 
24.5, and 25.5 are iteratively smoothed with a boxcar median filter of size 5, 15, and 21 pixels respectively. Lower 
right panel: A true c olor RGB composite i mage using the U, B, I filters. Each channel has been corrected for Galactic 
extinction following ISchlegel et al.l dl998f) and co nverted to the AB photometric system. The RGB composite was 
created by implementing the lLupton et al] (120041 ') algorithm. 
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Figure Bl. ES0185-13 A faint structure is visible to the North-East along the tail in all optical filters. This feature is not seen in the 
NIR, neither in H (not shown) or K. This could be tidal debris expelled from the main body of the galaxy, or a low luminosity dwarf 
which has passed through the galaxy. The average surface brightness in this feature is nu ~ 23.6, (j,b ~ 29.2, /iy ~ 28.9, fj-n ~ 27.7, 
and [ij ~ 27.8 mag arcsec -2 . The contrast and the intensity cuts have been selected to enhance the faint features but are identical in all 
filters to facilitate direct comparison. 




Figure B2. UM448 A faint structure is visible in the optical, extending in a tidal arm to the North-West. Though not as clear, parts 
of this feature are also visible in the NIR, North— West and South-East of the main body of the galaxy. The average surface brightness 
measured over the tidal arm is /ig ~ 26.5, Hv ~ 25.8, /ir ~ 25.5, and /ij ~ 24.5 mag arcsec -2 . The contrast and the intensity cuts are 
identical in all filters to facilitate direct comparison. 
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Figure B3. Extreme contours. Numbers in red/black indicate the isophotal level of the red/black contour in mag arcsec 
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